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ABSTRACT 

We present custom-processed ultraviolet, optical, and near-infrared photometry for the RESOLVE 
(REsolved Spectroscopy of a Local VolumE) survey, a volume-limited census of stellar, gas, and 
dynamical mass within two subvolumes of the nearby universe (RESOLVE-A and RESOLVE-B). 
RESOLVE is complete down to baryonic mass ~ iq9 i- 9.3 probing the upper end of the dwarf 
galaxy regime. In contrast to standard pipeline photometry (e.g., SDSS), our photometry uses optimal 
background subtraction, avoids suppressing color gradients, and employs multiple flux extrapolation 
routines to estimate systematic errors. With these improvements, we measure brighter magnitudes, 
larger radii, bluer colors, and a real increase in scatter around the red sequence. Combining stellar 
mass estimates based on our optimized photometry with the nearly complete HI mass census for 
RESOLVE-A, we create new z=0 volume-limited calibrations of the photometric gas fractions (PGF) 
technique, which predicts gas-to-stellar mass ratios (G/S) from galaxy colors and optional additional 
parameters. We analyze G/S-color residuals vs. potential third parameters, finding that axial ratio is 
the best independent and physically meaningful third parameter. We define a “modified color” from 
planar fits to G/S as a function of both color and axial ratio. In the complete galaxy population, 
upper limits on G/S bias linear and planar fits. We therefore model the entire PGF probability density 
field, enabling iterative statistical modeling of upper limits and prediction of full G/S probability 
distributions for individual galaxies. These distributions have two-component structure in the red 
color regime. Finally, we use the RESOLVE-B 2Icm census to test several PGF calibrations, finding 
that most systematically under- or overestimate gas masses, but the full probability density method 
performs well. 

Subject headings: galaxies: ISM — galaxies: photometry — surveys 


1. INTRODUCTION 

As imaging surveys provide ever more sky-coverage and 
greater depth, we are producing larger galaxy data sets 
probing to lower masses. Photometry from these imaging 
surveys allows estimation of stellar masses for galaxies, 
which only provides a partial view of galaxy mass with¬ 
out any cold gas data. The cold neutral gas mass probed 
by 21cm atomic hydrogen (HI) observations is generally 
the most abundant form of cold , observable gas in galax¬ 
ies in the nearby universe le.g.. lObreschkow &: Rawlingsl 
l2009f) . HI observations however can be time consuming, 
especially for galaxies with low absolute gas mass. 

Galaxies with low gas content can be of extremely 
different types: gas-poor galaxies of all stellar masses 
and gas-rich galaxies with low stellar masses. With ex¬ 
isting flux-limited surveys such as th e ALFALFA 21cm 
blind HI survey (|Havnes et al.l l2011h , we cannot mea¬ 
sure the gas masses for these two populations beyond 
our nearest neighbors. Fractional gas-mass limited sur- 
veys, such as the GA LEX Arecibo SDSS Survey (GASS; 
iGatinella et al.l[2fiTfll l and the Nearby Field Galaxy Sur- 
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vey (NFGS; iWei et aP l201d : iKannaopan et al.l I2013L 
hereafter K13), allow us to examine galaxy gas content 
for a wider range of galaxy types. Both of these data sets 
are representative of the galaxy population in that they 
sample all types of galaxies within their respective selec¬ 
tion criteria. Neither of these two samples, though, is a 
fair representation of the statistical distribution of galax¬ 
ies in the nearby universe. In contrast the RESOLVE 
(REsolved Spectroscopy of a Local VolumE) survey is a 
complete volume-limited data set that contains all galax¬ 
ies above a “cold baryonic” (stellar -|- cold gas) mass limit 
of -I091-9-3 Mq (in two separate subvolumes Kannap- 
pan et al. in prep.). The RESOLVE HI mass census is 
also fractional mass limited (Stark et al. in prep.). 

Already obtaining an HI mass census for the RE¬ 
SOLVE survey (^1550 galaxies) has required several 
hundreds of hours on radio telescopes. To obtain gas 
masses for larger galaxy data sets, we must develop ac¬ 
curate gas mass predictors. One particular use of such 
estimators is to obtain galaxy cold baryonic masses, 
which are the optimal indicator of dynamical mass for 
gas-r ich galaxies (e.g . , the baryonic Tully-Fisher rela¬ 
tion, iMcGaugh et al.l [2000f) . For higher mass galaxies 
the baryonic component is dominated by the stars. For 
lower mass galaxies, particularly below the gas-richness 
threshold mass at ^10®'^ Mq in stellar mass, galaxies 
can have as much cold gas as stars, or even be domi¬ 
nated by their cold gas mass (K13). It is important to 
characterize galaxy mass, especially for dwarf galaxies. 
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by cold baryonic mass (stars + cold gas) rather than stel¬ 
lar mass alone. For large imaging surveys, such charac¬ 
terization will be impossible without the aid of accurate 
gas mass predictors calibrated on existing galaxy surveys 
with complete HI data. 

One such predictor is the photometric gas fractions 
“PGF” technique, which allows us to estimate galaxy 
cold gas mass primari ly using colo r . The PGF technique 
was first presented in iKannapoanl (|2004ll as an observed 
relation between log g as-to-stellar-mass ratio or G/S and 
u — K color (see also iKannapoan fc Weil[2?)0^ . The re¬ 
lation between log(G/S) and color is surprisingly tight; 
cr ^ 0.37 dex. This early work on the PGF technique 
used a sample that cross-matched between a flux-limited 
parent sample from imaging surveys and a heterogeneous 
collectio n of available HI de tecti ons from the Hyper Leda 
catalog (jPaturel et al.ll200'^ . In iZhang et al.l (|2009[ 1 . the 
authors used a similarly selected sample and find smaller 
scatter tr ^ 0.3 dex using g — r color and including z-band 
surface brightness as a third parameter in the fit. 

More recently, the GASS team has explored the PGF 
technique using NUV—r color combine d with stellar mass 
surface density ((Gatinella et al.l 120101) to create a “gas 
fraction plane,” finding a = 0.315 dex. GASS is a stellar 
mass limited sample that is representative of high mass 
galaxies and has measured HI masses or upper limits 
down to a fixed fractional gas mass of 1-5% of the stellar 
mass. Their PGF calibration, however, does not accu¬ 
rately recove r the HI masses f or the blues t, most gas-rich 
galaxi es. In IGatinella et al.l (|2012ll and ICatinella et al.l 
(|2013ll . the authors provide updated calibrations exclud¬ 
ing galaxies with NUV—r > 4.5, which yield smaller 
residuals for gas-rich galaxies and smaller scatter over¬ 
all a = 0.29 dex. To co mbat the residuals for gas-rich 
galaxies, iLi et al.l (|2012t) use the GASS sample to pro¬ 
duced a calibration from a combination of NUV—r color, 
stellar mass, stellar mass surface density, and g — r color 
gradient. Their PGF calibration more accurately pre¬ 
dicts log(G/S) for gas-rich galaxies from the flux-limited 
ALFALFA survey with cr = 0.29 dex. The use of mul¬ 
tiple variables covariant with log(G/S) and each other, 
however, prevents meaningful physical interpretation and 
artihcially reduces scatter. 

The ALFALFA blind 21cm s urvey has also been used 
to derive a PGF calibration bv iHuang et al.l (|2012h . who 
use S/N > 6.5 reliable detections (code 1) and lower 
S/N detections with relia ble optical counterp arts (code 
2) from the a.40 catalog (|Havnes et al.l[2(Hltl . The cal¬ 
ibration is based on NUV—r color and stellar mass sur¬ 
face density. Since the ALFALFA survey is flux-limited, 
the calibration sample is biased towards gas-rich objects 
and produces an offset towards higher gas fr actions when 
comp ared to the GASS PGF calibrations (|Hnang et al.l 

[Ma). 

Lastly, K13 provide s a PGF calibration for the Nearby 
Field Galaxy Survey (I Jansen et al.ll2000tl . a H-band se¬ 
lected, representative galaxy survey that contains either 
HI detections or strong upper limits for all galaxies. The 
PGF calibration uses only u — J color and has scatter 
of a ~ 0.34 dex. While the scatter measured is higher 
than in other works, we note that the calibration relies 
on color only and includes low-mass galaxies, which have 
larger intrinsic uncertainties on their stellar mass esti¬ 
mates, while GASS is limited to high stellar mass galax¬ 


ies. K13 also shows the effect of adding molecular gas 
for a subsample of the NFGS galaxies, Hnding that for 
large spiral galaxies with low values of log(G/S) the cal¬ 
ibration is tightened when combining the molecular and 
atomic hydrogen mass as the galaxy cold gas mass. 

The interpretation of the tight relation between color 
an d 1 ok(G/S) has bee n discussed in a few of these works. 
In iKannappanl (1200411 the correlation between log(G/S) 
and u — K color is linked to the correlation between ap¬ 
parent u-band magnitude and apparent HI magnitude. 
This correlation is understood as the common link be¬ 
tween the two quantities and the amount of recent star 
formation within the galaxy. 

Another interpre tation of the PGF relation comes from 
I Zhang et al.l (1200911 . who claim the PGF calibration is a 
manifestation of the Kennicutt-Schmidt relationship be¬ 
tween the surface densities of star formation and of cold 
gas, which has been calibrated on the short s tar forma¬ 
tion t imescales probed by Ha (|Schmidtlll963l iKennicuttI 
II998I1 . The results of K13, however, show that the u — J 
color of a galaxy can be interpreted through stellar popu¬ 
lation modeling as the fractional stellar mass growth rate, 
defined as the mass of stars formed in the last Gyr di¬ 
vided by the pre-existing stellar mass. Thus u — J color 
probes timescales much longer than those probed by Ha. 
In this light, the current galaxy gas reservoir is related to 
the galaxy’s past growth rate over long timescales, and 
blue low-mass galaxies, which typically have high gas-to- 
stellar mass ratios (sometimes as much as 10), have been 
growing at rates inconsistent with closed box models and 
requiring ongoing cosmic accretion (K13). The authors 
argue that it is the long-term physics of accretion, rather 
than the short-term physics of the Kennicutt-Schmidt 
relation, that underlies the PGF technique. 

In this work, we provide new z=0 PGF calibrations 
using the A-semester of the volume-limited RESOLVE 
survey (RESOLVE-A). This data set offers several key 
advantages over the previous calibrations discussed here. 
First, we use newly reprocessed photometry, presented 
here, from several imaging surveys. Superior photometry 
and well understood systematic errors allows us to esti¬ 
mate reliable stellar masses through SED fitting. Second, 
we have an almost complete (78%) HI data set for galax¬ 
ies with detections or strong upper limits (defined here as 
IAMhi < 0.05Mstar), and we are able to incorporate the 
remaining 22% that are confused or have weak upper lim¬ 
its through statistical modeling using survival analysis. 
Third, our data set is limited on absolute r-band mag¬ 
nitude, which most closely corresponds to baryonic mass 
(K13), and the survey is complete to ^ 1 0^'^ Mq, 

well into the gas-dominated regime (see 112.11) . Lastly, 
because we use a volume-limited data set, we correctly 
represent the number density of galaxies in the local uni¬ 
verse in color and log(G/S) parameter space. 

This paper is organized as follows. First we describe 
the RESOLVE survey and its two subvolumes in 1J21 Next 
we detail the reprocessed photometry, stellar mass es¬ 
timation, and HI data in ^ We then describe color- 
limited PGF calibrations using linear fits in IjH and ex¬ 
amine correlations between log(G/S) residuals and pho¬ 
tometric parameters to obtain tighter calibrations in 1|51 
The linear hts are limited by their inability to predict gas 
masses for red galaxies, for which the correlation between 
color and log(G/S) breaks down, as well as by the fact 
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that we cannot simply include galaxies with weak upper 
limits. To properly predict gas masses for all galaxies, 
we describe in ^ a new PGF calibration using a 2D 
model to fit to a density field, yieldiirg log(G/S) prob¬ 
ability distributions for galaxies of all colors. In ^ we 
test the new calibrations on the RESOLVE-B data set 
and we compare with previous calibrations from the lit¬ 
erature, finding that our new calibrations are not biased 
for a z=0, volume-limited survey. Lastly we summarize 
our main conclusions in ^ 

2. DATA SETS 

For this work, we use the RESOLVE survey 
(|KannaDDan fc V^l2008t Kannappan et al. in prep.), a 
volume-limited mass census, to create and test new PGF 
calibrations. The RESOLVE survey is ideal for calibrat¬ 
ing gas mass estimators because it has a complete galaxy 
census with nearly complete HI data down to fixed frac¬ 
tional mass limits. 

RESOLVE is an equatorial survey covering two 
semesters (RESOLVE-A and RESOLVE-B) shown in 
Figure [T^. The RESOLVE survey is located within 
the SDSS footprint and makes use of the SDSS red- 
shift survey to build up survey membership with com¬ 
pleteness down to Mr,petro = —17.23, the absolute r- 
band magnitude corresponding to the SDSS survey limit 
of iQ-r^petro = 17.77 at the outer RESOLVE cz limit, 
7000 km s“^. We also include additional redshifts from 
various ar chival sources: t he Updated Zwicky Gata- 
logfUZC, iFalco et "nillQOOfl Hvoer Leda (iPaturel et al.l 
[200 . 6 dF (l.lones et al.ll2009ll . 2dF (IColless et al.ll2001ll 
GAMA ([Driver et al.l 1201 If) . ALFALFA ([Havnes et all 
M), and RESOLVE observations (Kannappan et al. 
in prep.). These extra redshifts provide greater com¬ 
pleteness to the RESOLVE data set, as detailed in a 
companion paper on the baryonic mass function and its 
dependence on environment (Eckert et al. in prep.) and 
in the RESOLVE survey design paper (Kannappan et 
al. in prep.). For both RESOLVE-A and RESOLVE-B 
we have custom reprocessed photometry providing to¬ 
tal magnitudes and systematic errors for GALEX NUV 
(plus new Swift UVOT imaging for nineteen galaxies), 
SDSS ugriz, UKIDSS YHK, and 2MASS JHK bands 
as available (described in 1 13.11) . 

To define survey membership, we use the redshift of 
the group to which each galaxy belongs. Group find¬ 
ing i s performed us i ng th e Friends-of-Friends algorithm 
from iBerlind et gdl (|2006[ ) with on sky and line of sight 
lin king lengths of 0.0 7 and 1.1 respectively as suggested 
by [Duarte &: MamonI l|2014[) and also justified in Eckert 
et al. (in prep.). As can be seen in Figure [1^, galax¬ 
ies with redshifts nominally outside the volume may be 
grouped with galaxies inside the volume, while occasion¬ 
ally galaxies with nominal redshifts inside the volume 
may be removed as they belong to a group outside the 
volume. 

The gas data used in this paper come f rom the RE¬ 
SOLVE HI census, which is described in 93.31 and will 
be published in Stark et al. (in prep.). RESOLVE HI 
observations build o n the ALFALFA blind 21cm survey 
([Havnes et al.l[2MT[) . which covers the entire RESOLVE- 
A footprint and partially covers the RESOLVE-B foot¬ 
print. New pointed observations with the GBT and 
Arecibo telescopes follow up on marginal detections. 


sources with weak upper limits, or sources with no HI 
data. 

2.1. RESOLVE-A 

The RESOLVE-A data set shown in Figure [T]a, 
occupies a volume of ^38,400 Mpc^ defined by: 
131.25° < RA < 236.25°, 0° < Dec < 5°, and 

4500 km s“^ < cz < 7000 km s“^. The data set’s r- 
band absolute magnitude distribution is shown in the 
orange solid line histogram in Figure [Ud. RESOLVE-A 
is complete down to < — 17.33 using the repro¬ 

cessed photometry described in 93.11 A magnitude of 
Mr,tot ~ —17.33 roughly corresponds to Mtary ^ 10®'^ 
Mq (K13). To determine the baryonic mass complete¬ 
ness limit, we consider the scatter in baryonic mass-to- 
light ratio, which can be at least as high as 3 resulting in a 
baryonic mass limit of 10®'^ Mq. The RESOLVE-A sur¬ 
vey contains 955 galaxies brighter than Mr,tot = —17.33. 
Of these 955 galaxies, ^12% were added from redshift 
surveys besides the SDSS main redshift survey. The data 
set resulting from the SDSS main redshift survey alone 
(RESOLVE-Aorig) is shown as the red dot-dashed line 
histogram in Figure [B d. The RESOLVE-A region is 78% 
complete in HI when counting successful, unconfused HI 
detections and strong upper limits resulting in I.AMhi 
< 0.05Mstar- We use the RESOLVE-A data set to de¬ 
termine our PGF calibrations, accounting for missing HI 
data with an iterative Monte Carlo technique akin to 
survival analysis (see 21and 

2.2. RESOLVE-B 

The RESOLVE-B data set is located in the 
SDSS Stripe 82 equatorial region, and it occu¬ 
pies a smaller volume of ~13,700 Mpc^ defined by: 
22h < RA < 3h, -1.25° < Dec < 1.25°, and 
4500 km s“^ < cz < 7000 km s“^. In Figure [Bd the 
absolute r-band magnitude distribution is shown for 
RESOLVE-B galaxies coming from the SDSS main red¬ 
shift survey as the black hashed histogram (RESOLVE- 
Borig), as well as for the full RESOLVE-B data set 
(grey filled histogram), which includes redshifts from 
the sources mentioned in and extra SDSS redshift 
observations over the Stripe 82 footprint. The data 
set is complete in r-band absolute magnitude down 
to Mr,tot — —17.0, slightly deeper than RESOLVE- 
A implying completeness to Mbary ~ 10®'^ Mq. The 
RESOLVE-B survey contains 487 galaxies to this limit, 
^25% of which have been added by redshift surveys be¬ 
sides the SDSS main redshift survey. We have recovered 
more galaxies in RESOLVE-B than in RESOLVE-A due 
to the extra spectroscopic passes done by the SDSS that 
are not part of the main SDSS redshift survey. The 
RESOLVE-B region is ~75% complete in HI data for 
good HI detections and strong upper limits. We use the 
RESOLVE-B data set to test our new PGF calibrations 
and compare with other calibrations from the literature 
(see [21). 

3. DATA 

For this work we need consistent and well calibrated 
photometry, stellar masses, and HI masses down to fixed 
fractional mass limits. We present our methods for re¬ 
processing UV, optical, and IR photometry for the RE¬ 
SOLVE survey in 93.11 We then describe our stellar mass 
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Figure 1. RA-cz and r-band absolute magnitude distributions of RESOLVE-A and -B semesters, a) The regions outlined in black show 
the RA and cz limits of the two RESOLVE subvolumes; both have been collapsed in Dec, ranging 0-5° for RESOLVE-A and —1.25 to 
+1.25° for RESOLVE-B. The black dots show RESOLVE galaxies, identified as members of RESOLVE because their group redshift falls 
within the limits of the survey. Orange points show galaxies belonging to groups outside of the RESOLVE cz limits 4500-7000 km s“^. 
b) The black hash filled histogram and red dot-dashed outlined histograms show the original absolute r-band magnitude distributions for 
RESOLVE-B and RESOLVE-A respectively. Both distributions fall off rapidly below Mr,tot < —17.33. The grey shaded histogram and 
the orange solid outlined histogram show the full absolute r-band magnitude distributions for RESOLVE-B and RESOLVE-A after redshift 
completion efforts described in m The RESOLVE-A region is still complete only to Mr,tot = —17.33, however we are able to move the 
RESOLVE-B completeness limit down to —17.0. 


estimation through SED modeling in 93.21 Lastly we de¬ 
scribe the various s ourc es of HI data and the measure¬ 
ment of HI masses 93.31 


3.1. Photometric Data 

We have reprocessed photometric data for the RE¬ 
SOLVE survey from the UV to near IR to obtain con¬ 
sistent, well-determined total magnitudes, and we use 
two to three methods of flux extrapolation per band 
to characterize systematic errors on the total magni¬ 
tudes of the galaxies. We have also run the same 
pipeline on the larger volume-limited ECO (Environmen¬ 
tal context) catalog (Moffett et al., submitted), which 
surrounds the RESOLVE-A subvolume. We us e opti¬ 
cal ugriz data fro m SPSS (lAih ara et aH l2011h . NIR 
JHK from 2M ASS (iSkrutskie et al.ll2006f) and/or YHK 
from UKIDSS (iHamblv et al.l 20081). and NUV from the 
GALEX mission (jMorrissev et alJ l2007f l. Our NUV 
data are mostly MIS depth due to prioritization of the 
RESOLVE-A footprint late in the GALEX mission (af¬ 
ter the FUV detector failed), while RESOLVE-B (Stripe 
82) already had deep coverage in both the NUV and 
EUV for other programs. The SDSS optical imaging in 
the RESOLVE-B footprint is extra deep due to repeated 
imaging with typic ally 20 frames per location on the sky 
(|Annis et al.ll20lll . With our improved photometry and 
realistic error measurements, we are able to measure reli¬ 
able colors and perform accurate stellar mass estimation 
via SED modeling. 

Our reprocessed photometry improves over SDSS 
pipeline photometry in several key ways. First, we use 
ima ges with improved sk y subtraction coming from ei- 
ther lBlanton et akl l)2011f) for SDSS or our own additional 
sky subtraction for 2MASS and UKIDSS. Second, we use 
the sum of the high S/N gri images to define the ellipti¬ 


cal apertures, allowing us to determine the PA and axial 
ratio of the outer disk if present. Third, we apply these 
same elliptical apertures to all bands which allows us to 
measure magnitudes for galaxies that may not have been 
detected by the original automated survey pipeline in 
certain bands, especially low surface brightness galaxies 
in 2MASS, UKIDSS, and GALEX. Lastly, we use two 
to three non-parametric methods of total magnitude ex¬ 
trapolation, measuring the light from each band inde¬ 
pendently (see Figure [Hand 93.1.2|) . This last point al¬ 
lows for color gradients within galaxies, as opposed to 
the model magnitudes provided by SDSS (more details 
in 93.1.31) , and allows us to measure systematic errors on 
magnitudes. 

We provide a comparison of the magnitudes, colors, 
and radii with photometry from the DR7 catalog of SDSS 
in Figure [3] and 93.1.31 Briefly summarizing, we find 
that the newly reprocessed photometry yields brighter 
magnitudes and larger effective radii. The colors tend 
to be bluer for large objects, which we believe to be a 
consequence of both the improved sky subtraction from 
[Blanton et al.l (120111) and the fact that we allow color gra¬ 
dients in magnitude estimation. The newly reprocessed 
photometry does not create a tight red sequence on the 
color-magnitude relation as seen in the DR7 photometry, 
however we argue that the tight red seque nce may be a 
consequence of these two issues in 93.1.31 We also dis¬ 
cuss an independent validation of our methods with the 
NFGS survey (shown in Figure 2a of K13) in 93.1.31 

In addition to the reprocessed photometry, this pa¬ 
per also presents new UV observations of 19 galaxies 
using the Swift Ult raviolet/ Optical Telescope (U VOT, 
iRoming et al.ll2005l . see also iGehrels et al.l I2004D . We 
use imaging from the uvm2 filter, which has a com¬ 
parable central wavelength but narrower width than 
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the GALEX NUV filter Isee iPoole et alJl2008f) . Com¬ 
pared to GALEX, the pointing restrictions for UVOT 
are much less stringent, allowing us to obtain obser¬ 
vations during Swift team fill-in time for RESOLVE-B 
galaxies that were not observed by GALEX or had only 
AIS depth (^150s) coverage. Nineteen galaxies were 
observed for more than 1 ks, the minimum exposure 
for useful photometry. Images were processed following 
iHoversten et al.l (1201 If) . Each galaxy was manually in¬ 
spected to make sure that the surface brightness in the 
uvm2 band was low enough that the resulting photom- 
etry errors due t o coincidence loss were below 1% (see 
Poole et ^ 120081 iBreeveld et al.l [2M^ IHoversten et al.l 
2011fl . We apply similar photometric processing to the 
Swift data as to the archival data including matched el¬ 
liptical apertures and multiple extrapolation techniques. 

3.1.1. Custom Processed Data 

We start the photometric reprocessing by download¬ 
ing the data from each respective website, performing 
background subtraction and coaddition when necessary, 
and cropping a region around the galaxy 9 times the 
Petrosian 90% light radius as reported by SDSS with 
a minimum crop size of 3x3 arcmin^. Because some 
galaxies are quite large, we rescale images to a fixed im¬ 
age size, causing the pixel scale to vary from galaxy to 
galaxy. Since the Stripe 82 region has been repeatedly 
observed in ugriz, for RESOLVE-B galaxies we coadd 
the many frames of data by inversely weighting by the 
variance of sky fluctuations using the IRAF task imcom- 
bine. Within the RESOLVE-A region there is typically 
one ugriz frame pe r region of sky, and we use SWARP 
(|Bertin et al.l 1200211 to stitch together adjacent frames 
when necessary, averaging together pixels where there 
is overlap between images. No additional background 
subtraction is done, as we are using SDSS DR8 im- 
ages with the optim ized sky background subtraction of 
iBlanton et all (DMHI . For 2MASS JHK and UKIDSS 
YJHK, we perform additional background subtraction 
by fitting and subtracting a 3rd order polynomial to a 
region of the galaxy frame where the galaxy and other 
objects are masked out. Coaddition is similar to the sin¬ 
gle frame SDSS process, using SWARP to stitch together 
2MASS and UKIDSS frames with a simple average to 
combine pixels in overlapping areas of the sky. Based on 
visual inspection of the UKIDSS data, we do not use the 
J band due to background subtraction and other issues 
that affect ^75% of the data. We have also examined 
the YHK images for each galaxy by eye to flag any cases 
with bad data. The GALEX NUV images do not require 
additional background subtraction, and these images are 
simply coadded using SWARP and weighted by exposure 
time. For the Swift uvm2 images we use Source Extrac¬ 
tor to identify and mask objects within the frame, then 
subtract off the median level of the non-masked areas as 
the sky-background. 

A significant number of galaxies (~16%) in RESOLVE 
have half-light diameters smaller than three times the 
typical r-band psf FWHM of ~1.4", warranting psf- 
matching across the optical bands and UKIDSS IR 
bands. For each galaxy, we use the SDSS provided ps- 
Field frame to reconstruct the psf for each band at the 
galaxy position on the frame. First we identify the band 
with the worst psf seeing (typically u or g). Next we find 


the Gaussian a value with which to convolve the psf of 
each given band to eliminate the difference between the 
psf of the worst band and that given band. This Gaus¬ 
sian a value is then used to create a Gaussian kernel 
that is convolved with the galaxy cropped image. Since 
we have changed the pixel scale of the frames of larger 
galaxies, we make sure to convert the a value into the 
correct pixel scale for that galaxy. For UKIDSS IR data 
a similar procedure is run, except that the psf for each 
frame is constructed from stars identified by Source Ex¬ 
tractor. If the value of the converted a value is less than 
one pixel, we do not perform the convolution. We do 
not psf-match the NUV, uvm2, or 2MASS JHK bands 
because their typical psfs are much larger than the SDSS 
(-5.5" for NUV, -2.5" for uvm2, and -2" for 2MASS). 
Thus aperture matched magnitude measurements for the 
UV and 2MASS IR will not be correct, especially for 
small galaxies. 

Masks ar e made from the r-band image using Source 
Extractor (|Bertin fc Arnouti 1199611 to detect stars and 
galaxies other than the target. Masks are checked by eye 
to ensure there is no under/over masking. In (13.1.31 we 
discuss the application of this photometric reprocessing 
pipeline for the ECO (Environmental COntext) catalog 
(Moffett et ah, submitted), for which we do not check 
each mask by eye. Instead we use a first iteration of 
the pipeline to check for discrepant total and aperture 
magnitudes as well as cases with no valid magnitudes at 
the end. For the most egregious outliers, we check the 
masks and edit by hand where necessary. 

To determine the parameters of the elliptical fit 
(namely the PA and ellipticity), we use an iterative pro¬ 
cedure involving two programs. First, Source Extractor 
is run on each galaxy’s cropped r-band frame to hnd an 
initial guess for the center, PA, and ellipticity, and 90% 
light radius. Second, we run the IRAF task ellipse on 
the gri coadded image using the Source Extractor quan¬ 
tities as inputs and allowing the PA and ellipticity, but 
not the center, to vary. The gri images have the highest 
signal-to-noise data, and by coadding these three bands, 
we provide the best image to feed to ellipse for deter¬ 
mining the PA and ellipticity in the the outer parts of 
the galaxy. The final PA and ellipticity are chosen using 
a median of the fits from the outer disk of the galaxy, 
where “outer disk” is defined between one and two times 
the 90% r-band light radius determined by Source Ex¬ 
tractor. Using the final PA, ellipticity and galaxy center, 
we then perform a fixed ellipse fit on the gri summed 
image to determine the set of annuli over which to mea¬ 
sure the galaxy surface brightness profiles for each band 
separately. 

These same annuli are then automatically applied to 
the GALEX, SDSS, 2MASS, and UKIDSS data. To 
match the NUV and YJHK images to the SDSS im¬ 
ages, we resample them to the pixel scale of the SDSS 
image. Imposing the same annuli over all bands allows 
us to measure the galaxy light out to its furthest extent 
(based on gri). For the IR bands, we are able to measure 
magnitudes for twice as many galaxies as the 2MASS cat¬ 
alog detects and for —1.15 times as many galaxies as the 
UKIDSS catalog (based on public DRSplus). 

Some galaxies are in very close pairs or embedded 
within a larger galaxy. To obtain better photometry for 
these galaxies, we have attempted to subtract off the 
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galaxy light from the interfering galaxy for each frame, 
when it seems possible to identify the light belonging to a 
specific galaxy. To flag such cases we search for whether 
a nearby galaxy is within 4 times the 50% r-band light 
radius of each RESOLVE galaxy. This selection returns 
27 systems. We inspect these 27 systems and remove 11 
which appear well separated from their neighbors. We 
also remove three closely paired systems (rsll58/rsll60, 
rsOlOO/rsOlOl, rs0196/rs0197), two that involve merging 
star forming spirals and one that contains two similarly 
sized elliptical galaxies, all three of which are so close as 
to make it impossible to disentangle the light from each 
galaxy. We find 13 systems that benefit from attempting 
to remove the light from a galaxy. 

Embedded Galaxies: There are five systems that are 
heavily embedded inside a much larger galaxy: rs0675, 
rs0749, rsl233, rsl227, and rs0072 (inside rs0673, rs0750, 
rsl232, rsl226, and rs0072 respectively). Another three 
systems are on the outskirts of a larger galaxy: rs0639, 
rsl089, and rf0090 (just outside of rs0642, rsl090, and 
rf0094 respectively). To obtain better photometry for 
these eight embedded galaxies, we first mask the small 
embedded galaxy and run ellipse on the larger galaxy. 
We then subtract off the model flux from the larger 
galaxy that is output from ellipse. We use the result¬ 
ing image that has the large galaxy subtracted out to 
run through the procedures described in this section, en¬ 
suring that any residuals from the model are masked out. 

Close Pairs: There are five close pair systems for 
which subtracting off the light of one or both of the 
members improves the magnitude estimates. These sys¬ 
tems are rs0267/rs0268, rs0397/rs0398, rs0851/rs0852, 
rf0015/rf0016, and rf0309/rf0310. To subtract off the 
light of each member we use the following steps: 

• First, we identify the galaxy with the simpler light 
profile, which we call galaxy-A. For each pair galaxy-A 
is rs0268, rs0397, rs0851, rf0016, and rf0310. 

• Second, we start with the image for the other galaxy, 
which we call galaxy-B. We mask galaxy-B and run el¬ 
lipse to fit the light profile of galaxy-A. 

• Third, we subtract off the galaxy-A model flux as 
provided by ellipse from the galaxy-B image. The result¬ 
ing image is used in the standard pipeline for galaxy-B 
(rs0267, rs0398, rs0852, rfOOlS, rf0309). 

• In most cases, we do not then subtract off the galaxy- 
B image for galaxy-A. We choose not to for a variety of 
reasons. For rs0267/rs0268, galaxy-B (rs0267) does not 
have a regular light profile making it difficult to sub¬ 
tract off. For rs0397/rs0398 and rf0309/rf0310 galaxy- 
B (rs0398, rf0310) is edge-on and easy to mask. For 
rs0851/rs0852, galaxy-B (rs0852) is much smaller and 
easier to mask out. 

• For the last pair rf0015/rf0016, we take the galaxy-B 
(rf0015) image with galaxy-A’s light subtracted off and 
mask any residuals from galaxy-A. We run ellipse to fit 
the light profile of galaxy-B. We then subtract off the 
model fit to galaxy-B provided by ellipse from the orig¬ 
inal galaxy-B image. The resulting image of galaxy-A 
(which is no longer in the center) is run through the 
pipeline with newly generated masks. 

We perform these procedures for all optical bands and 
UKIDSS images. For 2MASS and GALEX, we check first 
whether the subtraction is needed because the galaxy 
light may not extend far enough in these bands. 


3.1.2. Magnitude Extrapolation 

To extrapolate total magnitudes from the fixed ellipse 
fits for the optical bands, we use three methods: an ex¬ 
ponential (Sersic index n = 1) fit to the outer disk, a 
non-parametric Curve of Growth extrapolation, and an 
Outer Disk Color Correction based on the r-band. Fig- 
ure[2]shows schematics of all three methods in the g band 
for a RESOLVE-B galaxy. 

Outer Disk Eit: To compute the outer disk flux, we 
first define a fitting region where the annular flux is 1 to 
5 times the a of the sky noise. If the galaxy frame has 
been masked heavily (due to nearby stars or galaxies), 
we use a region 3 to 8 times the cr of the sky noise, 
and in extreme cases where a bright star or galaxy is 
on top of the galaxy, we use a region defined by 20-50 
times the a of the sky noise. Then we fit an exponential 
disk function to the fitting region of the galaxy surface 
brightness profile (between orange and red lines in left 
panel of Figure [2]) and sum the extrapolated flux from 
the inner edge of the last ellipse in the fitting region (red 
line) to an extremely large radius “Rex.” or 1000" past 
the inner edge of the last ellipse in the fitting region 
(red line). To compute the inner disk flux, we sum the 
raw, unmodeled flux interior to the inner edge of the 
last ellipse in the fitting region (red ellipse in left panels 
of Figure [2]). If pixels are masked in the raw data, we 
replace those values with the model output from ellipse. 
The exponential total magnitude equals the sum of the 
measured inner flux and the extrapolated outer disk flux. 
The typical dividing radius is near the 90% light radius 
in the r band. 

Curv e of Growth: The Curve of Growth method, fol¬ 
lowing iMunoz-Mateos et al.l (j2009f ) , computes the total 
magnitude using the derivative of the enclosed magnitude 
as a function of the radius. The enclosed magnitudes are 
calculated from the ellipse profile. A line is then fitted to 
the derivative of enclosed magnitude with respect to the 
radius vs. the total enclosed magnitude. The line is fit¬ 
ted over a new fitting region that extends farther out in 
the galaxy profile, to where changes in the total enclosed 
magnitude as a function or radius are small (between the 
light and dark green lines in the central panel of Figure 
HI). The y-intercept of the fitted line, where dm/dr = 0, 
is the total Curve of Growth magnitude. 

Outer Disk Color Correction: This method scales the 
outer disk r-band flux to determine the outer disk flux 
of an object in another band. First, we use either the 
Curve of Growth or the Outer Disk Fit r-band magni¬ 
tude to determine the radii containing 70% and 90% of 
the r-band light in the running total flux profile (light 
and dark blue ellipses show these respective radii in right 
most panel of Figure [2]). If the galaxy frame is heavily 
masked (more than 5% of the image pixels), we prefer 
the r-band exponential magnitude, otherwise the r-band 
Curve-of-Growth magnitude is used. We next measure 
the galaxy flux within the Rgo.r to Rgo^r annulus (Ao„t). 
If the S/N of the flux in this annulus is not greater than 
10, we decrease the inner radius of Aout by increments of 
5% down to 50% of the r-band light, stopping when we 
achieve S/N > 10. If the S/N is still less than 10 between 
Rgo.r and Rgo.r, we do not compute the galaxy magni¬ 
tude with this method. Otherwise, we calculate the flux 
ratio between a given band x and the r band within the 
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Figure 2. Illustration of the three methods described in 113.1.2l to extrapolate the total g-band magnitude for RESOLVE-B galaxy rf0218. 
In this case, all three estimates agree closely (17.55, 17.53, 17.55), yielding a small systematic error (0.021). Left column: To demonstrate 
the Outer Disk Fit method, we show in the top left panel annular g-band surface brightness vs. radius with the fitting region marked by 
the orange (inner) and red (outer) lines. The blue line shows the exponential disk fit to the data points. The bottom left panel illustrates 
how we compute the total magnitude as the sum of the raw galaxy flux inside the radius marked by the red line (same as radius in top 
panel marked by red line) and the extrapolated flux outside that radius. Middle column: To demonstrate the Curve of Growth method, 
we show in the top middle panel enclosed galaxy magnitude vs. radius with the fitting region marked by light green (inner) and dark green 
(outer) lines. The lower middle panel shows the enclosed galaxy magnitude vs. the derivative of enclosed magnitude with respect to radius 
for the points within the fitting region defined above. The blue line shows the fit and the y-intercept at dm/dr = 0 is the total galaxy 
magnitude. Right column: To demonstrate the Outer Disk Color Correction method, we show in the top right panel the annulus Aout 
defined by the r-band 70% and 90% light radii (light and dark blue lines respectively). Using the flux ratio in Aout, we fix the color of the 
galaxy to determine the g-band flux beyond Rgo,r, from the r-band flux in that region. 


annulus Aout, and we assume that this ratio continues 
out to infinity. From the r-band flux from Rgo.r to i?oo, 
and the flux ratio within Aout, we estimate the flux in 
band x from i? 9 o,r to Roo, then add this flux to the raw 
enclosed flux inside Rgo,r to get the final magnitude. 

Extrapolation of the NIR and UV magnitudes proceeds 
similarly to the optical extrapolation, but with a few 
subtleties. The Curve of Growth method is the preferred 
method for our NIR data due to the poor signal-noise for 
low surface brightness galaxies. Exponential fits are used 
to determine whether or not the Curve of Growth method 
works well. If the two hts disagree significantly or if the 
object’s magnitude is very faint, we look at the magni¬ 
tude calculated based on the i band (using the Outer 
Disk Color Correction method). If either the Curve of 
Growth or exponential matches the aperture magnitude, 
that is chosen. If neither method agrees, the Outer Disk 
Golor Gorrection is used and given a large systematic er¬ 
ror (>0.5 mag). For the UV data from both GALEX 
and Swift, we find that the Gurve of Growth method is 
the most reliable magnitude estimation method as the 
dumpiness of the UV and the possibility of XUV disks 
(extended U V emission outside th e typical optical extent 
of the galaxy; iThilker et al.l[200^ make exponential disk 
htting and hxing the outer disk color impractical. The 
Outer Disk Golor Gorrection method is also hampered 


by the mismatch in psf between the UV images (^5.5" 
and ^2.5" for the NUV and uvm2 respectively) and the 
psf of the convolved SDSS and UKIDSS images (^1.8"). 
If the Gurve of Growth method fails, though, we use the 
magnitude of the Outer Disk Color Correction method, 
with a systematic error > 0.06 applied. 

Errors for all bands are computed using not only the 
formal statistical error on the magnitude, but also the 
systematic error based on the difference in flux measured 
from the three methods. We apply a built in floor for 
the systematic error based on the overall distribution of 
systematic errors for the galaxy data set, such that none 
are lower than the original 25 percentile. 

In addition we compute half light and 90% light radii 
in the r band (i? 5 o,r and Rgo^r), as well as the r-band 
surface brightness within these radii {fir, 5 q and Rr, 9 o)- 
We also measure aperture magnitudes for all available 
bands within the r-band half light and 90% light radii, 
although the lack of psf correction for the 2MASS JHK 
and NUV and uvm2 bands compromises associated aper¬ 
ture matched colors. We also compute the g — r color 
gradient (hereafter Ag_j.), which is defined as the g — r 
color within the annulus between the half light and 75% 
r-band light radii minus the g — r color within the r-band 
half light radius. More positive colors indicate galaxies 
with bluer centers. 
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Throughout this work we use Milky Way foreground 
ex tinction correc ti ons d etermined from the dust maps 
of iSchle eel et ah (I1998D with the extinction curves of 
O’Donnelll (11994) fo r the optical and IR data, and of 
Cardelli et al.l (|19^ for the NUV and uvm2 data. For 
the NUV and uvm2 data we use the extinction correc¬ 
tion calculated at 2271 A and 2221 A, the effective wave¬ 
lengths of the NUV and uvm2 filter respectively. We 
note that using the more recently computed extinction 
coefhcients from ISchlafiv &: F inkbein ^ (|2011[1 . which use 
the extinction curve from iFitzpatricld (|1999ll . yields col¬ 
ors that are ~0.015 mag bluer in u — r (~0.04 bluer in 
NUV—r) and do not change the stellar mass estimates 
from 113.21 

Table [1] provides descriptions of the columns that are 
provided in a machine readable table with the photome¬ 
try for the RESOLVE survey. All galaxies processed are 
provided, including those in the buffer and fainter than 
the nominal RESOLVE limits. 

3.1.3. Comparison with Catalog Photometry 

We compare our newly reprocessed magnitudes, radii, 
and colors to the Petrosian and model photometry pro¬ 
vided in the SDSS DR7 catalog in Figure [3] SDSS cata¬ 
log Petrosian magnitudes are measured within a circular 
aperture of twice the Petrosian radius, defined as the ra¬ 
dius Rp where the ratio of the local surface brightness 
at Rn to the surface brightness within Rp is equal to 0.2 
(IBlanton et al.l[20Mll . The SDSS pipeline uses the Pet¬ 
rosian radius defined by the r band to compute Petrosian 
magnitudes for all other bands, thus yielding aperture- 
matched magnitudes and colors. The Petrosian system 
should pick up nearly total fluxes for disk (Sersic n = 1) 
galaxies, but is known to underestimate magnitudes for 
highe r Sersic n galaxies by ~0.2 mag (|Graham et al.l 
1200511 . The SDSS pipeline also computes model magni¬ 
tudes by fitting exponential (n = 1) and de Vaucouleurs 
(n = 4) models to the galaxy light profile, choosing the 
model of greater likelihood in the r band, and extrapolat¬ 
ing the profile to infinity. To measure model magnitudes 
for the ugiz bands, the SDSS pipeline scales the am¬ 
plitude of the r-band p rofile up or down to b est match 
the profile in that band (iStoughton et al.l[2002ll . Neither 
magnitude system is ideal as the Petrosian magnitudes 
do not measure the total galaxy light, while the model 
magnitudes are most sensitive to the inner profile of the 
galaxy and do not allow for color gradients within galax¬ 
ies. 

In Figure [3^ we compare our newly reprocessed r- 
band magnitudes with DR7 Petrosian catalog magni¬ 
tudes as a function of galaxy half light radius R^ow 
The reprocessed r-band magnitudes are overall brighter 
by ~0.13 mag than the DR7 Petrosian r-band magni¬ 
tudes. We find a similar, but slightly smaller, overall 
offset of ~0.1 mag between our newly reprocessed mag¬ 
nitudes and the DR7 model r-band magnitudes. The 
offset increases for the largest galaxies, as seen in the 
running median as a function of i? 5 o,T- (black dashed line. 
Figure [3]a,). Much of this trend can be attributed to our 
use of the improved sky background subtraction from 
IBlanton et al.l (|2011l) . which was not available for DR7. 
The blue solid line shows the expected median offset be¬ 
tween galaxy magnitudes using the new sky subtraction 
vs. the standard SDSS DR7 pipeline, as a function of 


Table 1 

RESOLVE Custom Photometry Catalog Description 


Column 

Description 

1 

RESOLVE ID 

2 

Right Ascension 

3 

Declination 

4 

cz 

5 

group cz 

6 

absolute SDSS r-band magnitude 

7 

apparent SDSS u-band magnitude 

8 

apparent SDSS w-band magnitude error 

9 

apparent SDSS p'-band magnitude 

10 

apparent SDSS p-band magnitude error 

11 

apparent SDSS r-band magnitude 

12 

apparent SDSS r-band magnitude error 

13 

apparent SDSS 2 -band magnitude 

14 

apparent SDSS i-band magnitude error 

15 

apparent SDSS z-band magnitude 

16 

apparent SDSS 2 :-band magnitude error 

17 

apparent GALEX NUV-band magnitude 

18 

apparent GALEX NUV-band magnitude error 

19 

apparent Swift uvm2-band magnitude 

20 

apparent Swift uvm2-band magnitude error 

21 

apparent 2MASS J-band magnitude 

22 

apparent 2MASS J-band magnitude error 

23 

apparent 2MASS iJ-band magnitude 

24 

apparent 2MASS iJ-band magnitude error 

25 

apparent 2MASS iC-band magnitude 

26 

apparent 2MASS iC-band magnitude error 

27 

apparent UKIDSS V-band magnitude 

28 

apparent UKIDSS V-band magnitude error 

29 

apparent UKIDSS Jf-band magnitude 

30 

apparent UKIDSS iJ-band magnitude error 

31 

apparent UKIDSS iC-band magnitude 

32 

apparent UKIDSS iC-band magnitude error 

33 

b/a axial ratio of outer disk 

34 

Rso^r half-light radius in r band 

35 

R 9 Q,r 90% light radius in r band 

36 

Ag—r g — r color gradient 

37 

{u — modeled u — r color 

38 

(u — i)'^ modeled u — i color 

39 

(u — J)^ modeled u — J color 

40 

{u — K)'^ modeled u — K color 

41 

{g — modeled g — r color 

42 

{g — i)'^ modeled g — i color 

43 

{g — J)’^ modeled g — J color 

44 

(g — K)'^ modeled g — K color 

45 

stellar mass 

46 

foreground extinction in u band 

47 

foreground extinction in g band 

48 

foreground extinction in r band 

49 

foreground extinction in i band 

50 

foreground extinction in z band 

51 

foreground extinction in NUV band 

52 

foreground extinction in uvm2 band 

53 

foreground extinction in Y band 

54 

foreground extinction in J band 

55 

foreground extinction in H band 

56 

foreground extinction in K band 


Note. — All magnitudes are newly measured from 
the raw images. Apparent magnitudes are provided with¬ 
out foreground extinction corrections. Foreground extinc¬ 
tion corrections used in this work are provided. Mod¬ 
eled colors designated by a superscript m are prod- 
ucts of the SED fitting routine from K13, described in 
EH and have foreground extinction corrections and k- 
corrections implicitly included. The datatable is provided 
at http://resolve.astro.unc.edu/data/resolve_phot_drl.txt 
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true galaxy Rsn.r (b ased on coefficients from Table 1 of 
iBlanton et al.lf2Mll only valid for i? 5 o,r > 5"). Our run¬ 
ning median matches very well with the expected trend. 
Note also that in this work, we do not use information 
from the inner profile of the galaxy to compute extrapo¬ 
lated total magnitudes, but rather extrapolate the light 
based on the outer profile of the galaxy. This difference 
may also contribute to the generally brighter magnitudes 
that we measure. 

In Figure [Sb we compare our newly remeasured i? 5 o,r 
values with DR7 catalog R^o.r values as a function of the 
newly remeasured i? 5 o,r- Since we measure greater flux 
per galaxy, we expect the half-light radii to be larger, 
and indeed we find that the new i? 5 o,r values are typi¬ 
cally ~49% larger than the SDSS Petrosian i? 5 o,r values 
and ^13% larger than the model i? 5 o,r values. The ratio 
between the new and catalog i? 5 o,r values becomes much 
greater a bove a remeasu r ed r of ~10", the value iden¬ 
tified by IBlanton et al.l (|2011li as the true galaxy half 
light radius above which the use of the new sky back¬ 
ground should significantly affect the measured galaxy 
flux and radius measurement. Another consideration af¬ 
fecting only the Petrosian radii is the fact that SDSS 
Petrosian apertures are circular whereas both our aper¬ 
tures and SDSS model apertures are elliptical. When 
we restrict the comparison of half light radii to galaxies 
with b/a > 0.85, our half light radii are only ^20% larger 
than the Petrosian radii, more in line with the 13% in- 
crease over the model radii. Figure 9 of lHall et al.l (|2012fl 
also shows the trend for Petrosian half light radii to have 
greater disagreement with remeasured half-light radii for 
more intrinsically edge-on galaxies. 

In Figure |4] we compare the newly reprocessed total 
u — r colors and the DR7 model u — r colors vs. Mr^tot- 
We find that the new u — r colors are overall ^0.18 mag 
bluer than the DR7 model u — r colors. To compute 
this offset, we measure the running medians of each color 
distribution as a function of M^yoi in 0.2 mag bins and 
subtract the two sets of median colors. We then de¬ 
termine the median of these median color differences, 
which is 0.18 mag. A large portion of the offset is 
due to the improved sky subtraction algorithm, but we 
also note the fact that our newly reprocessed photome¬ 
try allows for color gradients whereas SDSS model col¬ 
ors do not. Galaxy c olor gradients have been found in 
all galaxy types (e .g., Ide ■TonBlll99bl: l.lansen et al.l 1200(1 
ICibinel et ai.ll2013f)7 For example in the Nearby Field 
Galaxy Surve ^Jansen et al.l (|2000f) find that early and 
late types have typical B — R colors that are bluer by 
0.1-0.2 mag in their outer regions, while dwarf types dis¬ 
tribute evenly between blue and red color gradients. In 
this work we have explicitly allowed for color gradients by 
computing the total magnitudes in each band separately 
without the assumption of fixed profile shape built into 
the SDSS model magnitude algorithm. Even the Outer 
Disk Color Correction method fixes only the color outside 
the Aout annulus. 

A consequence of ignoring color gradients is that the 
red sequence defined by DR7 model u — r colors appears 
tighter than the red sequence defined by our newly re¬ 
processed u — r colors. To quantify the scatter, we fit 
a line to both sets of colors between the red sequence 
boundaries marked off by the red lines in Figure |4l and 
measure the rms from the fit. The red sequence defini¬ 


tion is shifted for the newly reprocessed u — r colors by 
0.18 mag to account for their overall bluer colors. We 
confirm the visual impression that the DR7 red sequence 
is tighter, finding that the SDSS model u — r red se¬ 
quence is tighter by ^16%. This tight red sequence seems 
to be an artifact of the SDSS model magnitude algorithm 
and should not be over-interpreted in measuring the star 
for mation historie s of re d sequence galaxies. We note 
that iSimard et al.l (|2011fl also report that using separate 
fits to compute g- and r-band magnitudes produces a 
more scattered red sequence than obtained when fixing 
fits in both bands to have the same half light radius, al¬ 
though these authors still choose fixed half light radius 
fits for convenience. 

In Figure[5]we compare independent photometric mea¬ 
surements for RESOLVE survey galaxies that overlap 
with the ECO (Environmental COntext) catalog (Mof¬ 
fett et ah, submitted), which is a larger volume-limited 
data set encompassing the RESOLVE-A subvolume. The 
ECO catalog has been reprocessed through the same 
pipeline, with the most significant difference in method¬ 
ology occurring at the masking step. Since ECO has ^10 
times the number of galaxies as RESOLVE, for ECO it 
was not feasible to check each mask by hand. The most 
egregious cases of over- or under-masking were deter¬ 
mined in a preliminary run of the photometry code on the 
catalog, by checking for extrapolated magnitudes that 
significantly disagreed with aperture magnitudes or cases 
where no magnitude was measured. The masks for these 
galaxies were then checked by eye to mitigate under- 
/over-masking. We compare the ECO and RESOLVE 
measurements for galaxies in the overlapping sub¬ 
volume in Figure [5^. We find no offset and only small 
differences of typically <0.2 mag between the two sets of 
magnitudes. Some of these differences may be attributed 
to the final magnitude chosen by the pipeline (Outer Disk 
Fit or Curve of Growth for the r band), which is based on 
the degree to which the frame is masked. FigurelSJ) shows 
the color-magnitude plots for both the full RESOLVE- 
A and RESOLVE-B (blue points) and ECO (orange-red 
contours) data sets, demonstrating that they are consis¬ 
tent. 

An independent validation of the methods used in this 
work is shown i n Figure 2a of K13 for the Nearby Field 
Galaxy Survey (|Jaiisen et al.ll2000h . All NFGS galaxies 
with available SDSS data were reprocessed through the 
same pipeline as described here and Figure 2a in K13 
shows that the reprocessed u — r colors are consistent 
with the expected Vega-AB offset ioi U — R colors mea¬ 
sured in iJansen et al.l l)2000f) over all angular sizes. The 
comparison between the total u — r color and the SDSS 
DR7 model u — r colors reveals an offset such that the 
new photometry yields ^0.2 mag bluer u — r colors than 
the DR7 model colors, similar to the offset that we mea¬ 
sure. 

3.2. Stellar Masses 

Stellar masses and k-corrected colors are calculated us¬ 
ing the spect ral energy distribution ( S EP) m odeling code 
described in iKannapoan fc Cawiseil (|2007[ 1. as modified 
by K13, which fits a grid of stellar population models 
to our newly reprocessed total NWugrizYJHK mag¬ 
nitudes plus new Swift uvm2 data for 19 galaxies). With 
photometric data from up to 10 bands, we are able to 
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Figure 3. Comparison of newly reprocessed photometry and SDSS DR7 photometry, a) Comparison of absolute Petrosian magnitudes 
y^r,petro and absolute total magnitudes Mr,tot from this work as a function of remeasured Our reprocessed photometry is brighter 

than the SDSS DR7 catalog Petrosian photometry; the running median is shown as a dashed black line. (We observe similar but slightly 
smaller offsets between our reprocessed photometry and the SDSS DR7 catalog model photometry.) The solid blue line shows the expected 
m agnitude difference us ing the new sky subtraction vs. the SDSS standard pipeline as a function of true galaxy based on the analysis 

of IBlanton et all H201in . which is in excellent agreement with our data, b) Comparison showing Rso.r from our reprocessed photometry 
divided by R 5 o,r from the DR7 photometry (Petrosian: black dots; model: pink squares), as a function of the new R 5 o,r in arcsec. The 
green solid line shows one-to-one correspondence, and the gree n dashed line marks 1 0^^ above which the new background subtraction should 
significantly affect the measured flux and radius of the galaxy lIBlanton et al.ir201in . Our newly reprocessed photometry yields ~13% larger 
half-light radii than SDSS model half-light radii, with larger increases for galaxies with R 5 o,r > 10'^ The increases over Petrosian radii 
are more extreme due to the assumption of circularity in the Petrosian algorithm. 


3.5 



-18 -20 -22 -24 -18 -20 -22 -24 

^r,tot ^r,tot 


Figure 4. Comparison oi u — r color vs. our Mr,tot for DR7 model colors (left) and our newly reprocessed colors (right). We determine 
the red sequence boundaries (red lines) for the DR7 data by eye, throwing out outliers that are too red. We then fit a line to the data 
points between the two red lines and compute the rms which is 0.12 for the DR7 model magnitudes (slope and width of red sequence shown 
by grey has marked region). To do the same for the newly reprocessed photometry, we first determine the overall color offset by computing 
the median color of each distribution in bins of 0.2 mag in Mr,tot. We then find the median of the differences of the medians in each bin, 
obtaining an overall shift of 0.18 mag. We shift down our red sequence boundaries by 0.18 mag for the newly reprocessed photometry. We 
then fit a line and compute the rms to be 0.14 mag (slope and width shown in green solid region with slope and rms of DR7 model colors 
shi fted d own by 0.18 mags and overplotted), which is larger by 0.02 mag or 15% than the rms computed for the DR7 colors. As argued 
in ^3.1.31 we believe the higher scatter in our red sequence to be more correct, as our photometry does not suppress color gradients. Also 
most of the few extremely red outliers (including one at u — rtot = 3.6 off the plot) are low surface brightness red galaxies whose shallow u 
band data are the hardest to measure. A few others are those embedded galaxies for which we subtracted off the light of the larger galaxy. 
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Figure 5. Comparison of photometric measurements for galaxies in both the RESOLVE survey and the ECO catalog, a) Difference 
in Mr,tot measured for the same galaxies in RESOLVE-A and ECO vs. the RESOLVE-A Mr,tot- The main difference between the two 
measurements arises in the masking step. Every RESOLVE-A galaxy mask is checked by eye, but for the much larger ECO data set 
we check only images for which we have identified large discrepancies between extrapolated and aperture magnitudes. The two sets of 
measurements agree well, with differences mostly <0.2 mag. Differences become larger for fainter galaxies, for which we expect larger 
uncertainties in extrapolation, b) Color-magnitude relations of the full ECO (contours) and RESOLVE (dots) data sets, show that they 
are consistent. While we do not use ECO in this work, the ECO catalog an extension of the RESOLVE survey, so it is useful to establish 
that its photometry is consistent with the rest of RESOLVE. 
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estimate robust stellar masses. We omit UKIDSS YHK 
values if the frames have been flagged by eye. We also 
omit UKIDSS HK and 2MASS JHK if the values are 
fainter than 18, 17.5 and 16, 15, 14.5 respectively. We 
also remove any NUV magnitudes fainter than 24, and 
we remove the u band magnitudes for four galaxies for 
which the u band data available from SDSS are essen¬ 
tially frames of noise. 

In this work use the second model grid from K13, 
which is a grid of composite stellar population models 
(CSPs) including an old simple stellar population (SSP) 
ranging in age from 2-12 Gyr and a young population 
either described by continuous star formation starting 
1015 Myr ago and turning off between 0 to 195 Myr ago 
or as a quenching burst with SSP age 360, 509, 641, 806, 
or 1015 Myr. The contribution from the young pop¬ 
ulation ranges from 1-94.1% of the stellar mass. The 
model g rid is built using the ste llar population mod- 
els from |Bruzual &: Chariot (|2003f) with a Chabrier IMF 
(|Chabried 1200311 . and four possible metallicities (Z = 
0.004, 0.008, 0.02, or 0.05). Eleven reddening values 
(rv ranges from 0-1.2) are applied to the young popula¬ 
tion fo llowing the dust extinction law given by ICalzettil 
(|2001ll . There is no physical or spatial model assumed for 
the dust, only an empirical determination of the amount 
of reddening and extinction based on the stellar popula¬ 
tion model grid fits to the galaxy SED. 

To determine a galaxy’s stellar mass, the stellar mass 
is computed for each CSP model in the grid and given 
a likelihood based on the value of the model fit to 
the data. Combining likelihoods over all models yields a 
stellar mass likelihood distribution for each galaxy. The 
median value of this stellar mass likelihood distribution 
is taken to be the nominal stellar mass of the galaxy. 

The SED modeling code also outputs the likelihood 
weighted colors for each galaxy, which are effectively 
“smoothed” by the model fits and implicitly k-corrected. 
We denote the use of these modeled colors with a super¬ 
script m (following the notation from K13 and Moffett et 
ah, submitted). The stellar population code also outputs 
de-extincted galaxy magnitudes, taking into account the 
internal extinction due to dust in the galaxy. These mag¬ 
nitudes cleanly divide the red and blue sequences in the 
color-stellar mass diagram as shown in K13 and Mof¬ 
fett at al., submitted. Here, however, we choose to use 
the modeled colors, which represent the actual rest-frame 
colors of the galaxies, for easier application of the PGF 
calibrations to other data sets. We show the color-stellar 
mass diagram for RESOLVE-A and RESOLVE-B using 
both u — r total colors measured from the raw repro¬ 
cessed photometry and {u — r)™ colors from the model 
hts in Figure El The SED modeled colors agree with 
the measured colors well within the expected k-correction 
values at these redshifts of up to 0.03 mag for the r band 
and 0.1 mag for the u band. 

Since some stellar mass estimation techniques have 
been shown to be biased as a function of inclination 
(|Maller et al.ll200^ . it is important to test whether our 
stellar masses may be biased as a function of axial ratio. 
We perform two tests. First, we select only galaxies with 
gas-to-stellar mass ratio > 0.1, implying significant gas 
and thus potentially dust, and we divide this subset into 
quartiles based on their axial ratio. A KS test reveals 
that the stellar masses of the upper and lower quartiles 
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Figure 6. Color vs. stellar mass for the entire RESOLVE survey. 
Black dots show total reprocessed magnitudes u — r and smaller 
red dots show the SED modeled colors (u — r)'^. While there are 
outliers in the reprocessed photometry as seen in Figure[4] the SED 
modeled colors use information over the entire NUYugrizYJHK 
SED and provide a cleaner color-mass plot. 

{b/a > 0.77 and b/a < 0.39) are consistent with being 
drawn from the same population {pnuii = 0.99). Second, 
we recompute our stellar masses applying the dust law 
to both the young and old populations (as opposed to 
just the young population as for our preferred mass esti¬ 
mates) . We find a tiny overall offset for late type galaxies 
of ~0.02 dex but no systematic trend between the two 
stellar mass calculations as a function of axial ratio. For 
early type galaxies we find a tiny differential systematic 
offset of 0.02 dex between the most elongated and round¬ 
est galaxies. Both of these tests suggest our stellar mass 
calculations are not biased by dust extinction. 

3.3. HI Masses 

The HI masses and upper limits for RE SOLVE come 
from the blind 21cm ALFALFA survey (iHavnes et al.l 
1201 If) and our own new observations with the GBT and 
Arecibo telescopes. It is important for creating a gas 
mass estimator to have complete HI data for the entire 
data set. Below we describe the observations taken to 
date, how we determine and handle confused sources, 
and the HI completeness of the RESOLVE data set. 

The ALFALFA survey has covered the entire 
RESOLVE-A region and the Dec 0° to 4-1.25° strip 
of RESOLVE-B, providing HI detections or upper lim¬ 
its (not necessarily strong) for 85% of RESOLVE. 
Data reduction and source extraction are described in 
IHavnes et al.l (j2011l) . At the nominal S/N limit of 6, the 
ALFALFA flux limit translates to a fixed HI mass sen¬ 
sitivity at RESOLVE distances of ^10® Mq. Since RE¬ 
SOLVE galaxies range from 10®-10^^'^ M 0 , this fixed 
sensitivity implies a large number of upper limits that 
are much weaker than our stated goal of IAYIhi < 
O.ObMstar- To increase the yield from the basic AL¬ 
FALFA data products. Stark et al. (in prep.) extract 
140 lower S/N detections and upper limits for RESOLVE 
galaxies within the ALFALFA grids. 

To further increase the useful HI data set, we 
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have acquired pointed observations with the GET and 
Arecibo telescopes obtaining HI data for 290 galaxies in 
RESOLVE-A and 337 galaxies in RESOLVE-B (Stark et 
al. in prep.). We target galaxies with either no HI mea¬ 
surements or weak upper limits from ALEALFA, aiming 
for detections with S/N ^ 10 or strong upper limits. In 
addition, we have cross -matched the RESOLV E catalog 
with the HI catalog of lSprineob et al.l (|2005ll to obtain 
thirteen more HI measurements. 

To check for consistency between our GET and Arecibo 
pointed observations, we have remeasured HI fluxes for 
^10 galaxies in RESOLVE and find consistency be¬ 
tween observations with the two telescopes within ^15- 
20% (Stark et al. in prep.). Gonsistency checks be- 
twe en ALFALFA an d Arec ibo pointed observations from 
the ISpringob et ^ (|2005ll catalog are documented in 
iHavnes et al.l (I20TI ) and HI flux measurements between 
the two catalogs are shown to be in agreement within 
~20%. 

HI masses and upper limits are calculated as described 
in Stark et al. (in prep.) and KI3. Confusion is deter¬ 
mined based on the source of the HI measurement, 4' 
for the smoothed resolution element of ALFALFA, 9' for 
the GET, and 3.5' for Arecibo pointed observations. De¬ 
confusion is performed following techniques described in 
Stark et al. (in prep) that improve on methods described 
in K13. The de-confused HI masses are provided in Stark 
et al. (in prep.). Neutral gas masses are calculated by 
multiplying the HI mass by 1.4 to account for helium, 
Mgos = \ AMhi- 

For this work, we apply the following set of criteria 
to determine reliable HI masses. We require detections 
to have S/N > 5. We use de-confused HI masses if the 
systematic error on the de-confused HI masses is < 25% 
of the de-confused HI mass. For limits, we require that 
the upper limit yields a gas mass < 0.05Ms4ar- 

Based on these criteria, we provide the statistics on 
HI completeness for this work for the two RESOLVE 
subvolumes. RESOLVE-A has a total of 955 galaxies, 
of which 637 have reliable HI detections with S/N > 5 
(34 of those detections are successfully deconfused ob¬ 
servations) and 107 have strong upper limits resulting in 
Mgos < O.OSMgtar- Thus 78% of the sample (744 galax¬ 
ies) have reliable HI data for defining PGF calibrations. 
In RESOLVE-B there are 487 galaxies, of which 294 have 
good detections (34 are successfully deconfused observa¬ 
tions) and of which 70 are strong upper limits, yielding 
75% of RESOLVE-B or 364 galaxies that have reliable 
HI data. 

In RESOLVE-A, we are still lacking adequate HI mea¬ 
surements for 211 galaxies. Of those 211 galaxies, 92 are 
weak upper limits yielding gas masses that range from 
0.052 - 4.01 X Mstar, 16 are low S/N detections, and 
103 have HI profiles where de-confusion is not possible. 
To examine whether these galaxies with inadequate HI 
measurements are biased, we plot color vs. log(Mstar) 
for RESOLVE-A in Figure 0 Galaxies for which we 
have weak upper limits (red dots) tend to be low-mass, 
red galaxies which are generally gas-poor and require the 
longest observing times for successful detection or strong 
enough upper limits. Galaxies for which we have low 
S/N detections (green dots) are low-mass blue objects or 
higher mass red objects. Galaxies for which de-confusion 
is impossible (blue dots) are scattered throughout color 



log M3,,, 


Figure 7. Color vs. stellar mass relationship for the RESOLVE-A 
data set. Black open circl es show galaxies with reliable HI detec¬ 
tions as described in 3 m] Red dots show galaxies with HI upper 
limits yielding Mgas > O.OSMstar, green dots show galaxies with 
low S/N detections (S/N < 5), and blue dots show galaxies with 
HI profiles that are impossible to de-confuse. The weak upper lim¬ 
its tend to be the low-mass red galaxies, which are generally gas 
poor and have the least absolute HI content of galaxies in RE¬ 
SOLVE. The low S/N galaxies are also low mass, but generally 
bluer. The confused galaxies are interspersed throughout color 
and stellar mass. 

and stellar mass. 

4. COLOR-LIMITED PGF CALIBRATIONS 

In this section we describe our method to provide z=0 
PGF calibrations via linear fits between log(G/S) and 
color. In Figure [SJ we show the relationship between 
log(G/S) and (m — J)"* color, which is clearly linear. 
However, for galaxies redder than (u — J)"* = 3.6 mag 
there is a breakdown in the correlation. While the cor¬ 
relation between log(G/S) and {u — J)™ color continues 
for some galaxies redder than 3.6 mag, we also see that 
the population of quenched galaxies with very low val¬ 
ues of log(G/S) becomes more important for these same 
red colors. In we describe a new calibration method 
using a 2D model ht to the probability density field of 
log(G/S) vs. color that allows us to model all galaxies. 

Here we are generating linear fits to predict values of 
log(G/S) from color, where the latter has much smaller 
errors and thus functions as a classical independent vari¬ 
able. Thus, especially given the likelihood of the in¬ 
trinsic scatter over and above the errors, to obtain the 
best p redictor we should minimize residuals in lo g:(G/S) 
alone (llsobe et al.l 119901 : iFeigelson fc Babul 1199211 . Due 
to the population of red galaxies with strong HI upper 
limits, we must exclude all galaxies redward of a ver- 
tical color cutoff, e.g ., (it — J)"* > 3.6 mag, similar to 
iCatinella et al.l (12012 ) and K13. Making a cut in color is 
appropriate for measuring the correct calibration to pre¬ 
dict log(G/S) from color as we want to preserve the scat- 
ter for the predicted quantity ([Kannapnan et aH I2002D 
rather than fitting to only the HI detections or making a 
cut in log(G/S). Excluding red galaxies, however, limits 
the validity of our PGF calibration to galaxies blueward 
of the red color cutoff. 

A set of such color-limited PGF calibrations for a va- 
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(U-J)"" 


Figure 8. Photometric gas fractions relation for u — J SED mod¬ 
eled color, hereafter (it — J)™'. For the photometric gas fractions 
relation gas refers to the atomic gas content, for which we in¬ 
clude a helium correction factor: NLgas = When fit¬ 

ting, we exclude galaxies redder than (it — J)"* = 3.6 mag, where 
the quenched galaxy population overwhelms the low log(G/S) end 
of the relationship between log(G/S) and color. We also exclude 
galaxies bluer than 2.0 mag as there are only a few data points 
that may skew the overall fit. 

riety of color combinations is summarized in Table [2j To 
create these, we use the 744 galaxies from the RESOLVE- 
A data set that have reliable HI detections or strong up¬ 
per limits. “Reliable” HI detections are considered to 
include non-confused detections with S/N > 5 as well as 
de-confused detections where the systematic error on the 
deco nfuse d gas mass is < 25% of the measured gas mass 
(see iSSl). We define strong upper limits to be those for 
which the gas mass is < 5% of We exclude galaxies 

redder than the red color cutoff of each color distribution 
listed in Table [2] (roughly where the upper limits start 
to dominate). We also trim points at the blue end where 
the density of points is low and outliers may affect the 
fit as indicated in Table [2] For {u — J)"* color, the blue 
color trim is 2.0 mag and the red color cutoff is 3.6 mag. 
Finally we perform an ordinary least squares forward fit 
to minimize the scatter in log(G/S), the quantity that 
we want to predict. We choose not to weight the fit by 
the measurement uncertainties in log(G/S), because they 
are correlated with the values of log(G/S) and color, so 
weighting by them would bias the fits towards galaxies 
with high gas content. The slope and offset in log(G/S) 
of these color-limited PGF calibrations are given in Table 
[2] along with the measured scatter in the relations and 
the blue color trim and red color cutoff valuesO 

We expect these fits to be useful for galaxies blueward 
of the blue color trim (as discussed in ^ but these cali- 

^ We note that the RESOLVE-A region, which we use for these 
linear fits, is less redshift complete than RESOLVE-B. We have 
performed empirical completeness corrections based on luminosity 
and surface brightness or color for the ECO catalog (Moffett et 
ah, submitted), which encompasses the RESOLVE-A subvolume. 
These empirical completeness corrections are based on the more 
complete RESOLVE-B subvolume. We find that weighting the lin¬ 
ear fits by these completeness corrections does not change the lin¬ 
ear fit parameters significantly and we do not use the completeness 
corrections in this work. 


brations do not allow us to predict gas masses for galaxies 
redder than the red color cutoff. We also note that color- 
limited linear calibrations, and all calibrations based on 
simple fits, are subject to bias without survival analy¬ 
sis to model galaxies that are confused, have weak up¬ 
per limits, or lack reliable HI detections. Routines to 
incorporate upper limits in a linear fit exist but rely on 
the assumption that the upper limits are dist ributed ran¬ 
doml y throughout the sample (as discussed in llsobe et al.l 
119861) . Such “random censoring” is not the case for the 
PGF calibration, as those galaxies with upper limits in 
the RESOLVE-A data set are primarily red galaxies with 
low gas-to-stellar mass content. Thus using such routines 
would not be statistically robust. In )j6]we improve on 
simple fitting by producing a 2D model of the log(G/S) 
vs. color probability density field. This fully probabilis¬ 
tic approach allows us to implement a version of survival 
analysis that reinserts the galaxies left out of the lin¬ 
ear fits and to predict log(G/S) probability distributions 
for individual galaxies, even those redder than the color 
cutoff. Before performing this analysis, however, we ana¬ 
lyze whether the residuals from these linear fits correlate 
with any other photometric parameters that may help 
produce tighter PGF relations. 

5. CORRELATIONS WITH 3RD PARAMETERS 

In this section, we seek a combination of color and 
other photometric parameters that may produce a tighter 
PGF calibration for more accurate gas mass estimation. 
To this end we use the RESOLVE-A data set to ana¬ 
lyze correlations between various photometric parame¬ 
ters and residuals from the color-limited PGF calibra¬ 
tions in T5.1I We then explore possible physical reasons 
for these residual correlations by examining their relation 
to galaxy morphology in f|5.2l Lastly we provide plane 
fits between color, axial ratio, and log(G/S) for tighter 
color-limited PGF calibrations in T5.3I 

5.1. Best 3rd Parameter for Gas Mass Estimation 

As potential third parameters, we examine the sur¬ 
face brightness within i? 5 o,r (/Tr.so): concentration index 
(Cr), photometric axial ratio {b/a), g — r color gradient 
(Ag_r), stellar mass and baryonic mass Mbary 

The parameter gr ,50 is found by determining i? 5 o,r from 
the ellipse profiles produced in the photometric repro¬ 
cessing. The r-band light within i? 5 o,r is then divided 
over a circular area defined by i? 5 o,r, giving a somewhat 
intrinsic surface brightness for each object, although we 
have not applied a correction for internal galaxy extinc¬ 
tion or optical depth since the goal is to provide a simple 
empirical recipe for gas estimation. Concentration in¬ 
dex is defin ed as Ran r /Rbn r follo wing IShimasakn et al.l 
(j2001|) and iStrateva et aLl (|2001h . The b/a measure¬ 
ment comes from the fixed ellipse fits in our photometric 
pipeline, ^g-r is defined as the g — r color within the an¬ 
nulus between the r-band 50% and 75% light radii minus 
the g — r color within the r-band half light radius, with 
more positive n umbers meaning th e gala xy has a bluer 
center following iKannappan et ahl (j2004f l. The stellar 
mass Mstoj. is the median stellar mass from the likeli¬ 
hood weighted SED models described in ^3.21 and the 
baryonic mass is Mstar + ^gas- 

Figure [9] shows log(G/S) residuals from the {u - J)^ 
color-limited PGF calibration plotted against these dif- 
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Figure 9. Measured minus predicted log(G/S) for RESOLVE-A from the PGF calibration for (u — J)'^. Residuals are plotted against a) 
r-band 50% surface brightness /i.r,50i b) concentration index Cr (Rgo/Re), c) axial ratio b/a, d) color gradient Ag—r^ e) stellar mass Mstar, 
and f) baryonic mass Mbary We use the Spearman Rank test to assess whether a correlation exists between the residuals in log(G/S) 
and these photometric parameters. The probability of no correlation and the strength of each correlation are reported for each panel. The 
parameters firSO, Ag—r, and b/a have significant correlations with the residuals in log(G/S). The correlation with Mstar is not physically 
meaningful but is caused by covariance as well as our selection on Mr,tot, a proxy for baryonic mass, which implies that the data set 
contains low stellar mass gas-rich objects but not low stellar mass gas-poor objects. We mark in panel e the stellar mass completeness limit 
for RESOLVE-A (M^ary = 9.0), below which the residuals are clearly biased towards higher log(G/S) measured than predicted. We also 
mark in panel f the baryonic mass completeness limit for RESOLVE-A (Mbary = 9.3). Below this limit we lack extremely gas-dominated 
residuals, which is expected due to the absolute magnitude limited sample, where the highest high baryonic mass-to-light ratio will fall 
below our Mr,tot = —17.33 mag limit. 
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Table 2 

Color-limited Photometric Gas Fraction Calibrations 


color 

(mag) 

slope 

(dex/mag) 

log(G/S) offset 
(dex) 

(7 

(dex) 

blue trim 
(mag) 

red cutoff 
(mag) 

N galaxies 

{u - r)"* 

-1.763 

2.725 

0.319 

1.0 

2.0 

552 

(u - i)^ 

-1.421 

2.510 

0.314 

1.0 

2.3 

571 

(u - J)"* 

-1.127 

3.337 

0.322 

2.0 

3.6 

560 

(u- K)”^ 

-1.059 

3.993 

0.331 

2.8 

4.4 

543 

(g - rT 

-3.488 

1.467 

0.302 

0.1 

0.6 

568 

(9 - ir 

-2.399 

1.546 

0.310 

0.2 

0.9 

557 

(9 - JT 

-1.582 

2.918 

0.332 

1.1 

2.2 

550 

(9 - 

-1.401 

3.744 

0.364 

2.0 

3.0 

501 
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ferent parameters. Only the data from galaxies used in 
the linear ht are shown, i.e., those having a color between 
the blue trim and red cutoff colors and reliable HI data 
or a strong HI upper limit. We have performed Spear¬ 
man Rank tests to assess whether there is a correlation 
between the residuals in log(G/S) and these third pa¬ 
rameters and quantify the strength of that correlation. 
Third parameters ^r, 50 j ^g-rj and all show 

some correlation with residuals in log(G/S), while Cr and 
Mbary show no significant correlation. Below we examine 
each significant correlation to determine the best param¬ 
eter for more accurate gas mass prediction. 

The correlation between residuals in log(G/S) and 
Mstar, shown in Figure [Qj:, is in the sense that low stel¬ 
lar mass objects tend to have higher measured log(G/S) 
than predicted, which is expected from covariant errors in 
log(G/S) and Furthermore, the bias becomes very 

evident for stellar mass galaxies <10® ° Mq (roughly the 
limiting stellar mass completeness limit) marked in panel 
e by a dashed line. This correlation between log(G/S) 
residuals and stellar mass is caused by our survey selec¬ 
tion o n Mr.tnt, which serves as a close proxy for baryonic 
mass (|KannaDDan fc Weill2008L K13) rather than stellar 
mass. For a volume-limited, baryonic mass limited data 
set, an abundance of gas-rich low stellar mass objects is 
built in. The correlation between log(G/S) residuals and 
stellar mass is thus not surprising and only reflects the 
r-band magnitude selection. In Figure IHf we see that the 
color-limited PGF calibration does not show evidence for 
a correlation between residuals in log(G/S) and baryonic 
mass. 

The remaining correlations between log(G/S) residu¬ 
als and photometric parameters that we investigate are 
Mr, 50, b/a, and Ag-r- These correlations may be re¬ 
lated to galaxy morphology and/or ev olut ionary state, 
which we explore in Figure [TO] and 115.21 To decide 
which of these third parameters is best for use in our 
PGF calibrations, we use the following criteria: I) the 
third parameter has a strong correlation with residuals 
in log(G/S), 2) the correlation between the third parame¬ 
ter and log(G/S) residuals is not covariant with the built 
in correlation between residuals in log(G/S) and stellar 
mass, and 3) the third parameter is a reliable, well de¬ 
termined quantity for all galaxies in the data set. 

First we examine Mr,so in Figure |9li. The correlation 
between residuals in log(G/S) and Mr ,50 is such that lower 
surface brightness galaxies have larger log(G/S) values 
than predicted. The Spearman Rank test shows that the 
strength of the correlation between residuals in log(G/S) 
and fir ,50 is high (R = 0.43) with small probability of 
no correlation (P = 3x10“®^). To test whether stellar 
mass is impacting the correlation between residuals in 
log(G/S) and Mr, so, we fit a line to the correlation be¬ 
tween residuals in log(G/S) as a function of stellar mass, 
and remove the stellar mass dependence from the resid¬ 
uals. We then run a Spearman Rank test on the cor¬ 
relation between the stellar mass corrected residuals in 
log(G/S) and Mr, so, finding that the strength decreases to 
R = 0.28 and the probability of no correlation is larger 
by several orders of magnitude (P = 9x10“^®) though 
still significant. The reliability of Mr,50 values is variable 
since Mr, so is evaluated within the half light radius. For 
the smallest galaxies, the quantity may not truly repre¬ 


sent the surface brightness within i? 5 o,f 

Next we examine the third parameter b/a shown in 
Figure[9t. The correlation between residuals in log(G/S) 
and b/a is such that more edge-on or disky galaxies have 
larger log(G/S) values than predicted. The strength of 
the correlation between residuals in log(G/S) and b/a 
is also high (R = —0.39) with small probability of no 
correlation (P = 3x10“®®). Running a Spearman Rank 
test on the correlation between the stellar mass corrected 
residuals in log(G/S) and b/a, we find that both the cor¬ 
relation strength and the probability of no correlation 
remain mostly the same (R = —0.37 and P = 5x10“®°) 
showing that the correlation is not affected by 
The b/a measurements are reliable, since they are evalu¬ 
ated over the outer disk of the galaxy. For the smallest 
galaxies there may be a tendency to measure rounder 
objects, but since b/a is evaluated in the outer disk, this 
issue should be minimized. 

Lastly we consider the third parameter Ag-r shown in 
Figure[9ji. The correlation between residuals in log(G/S) 
and Ag-r is such that galaxies with larger values of Ag-r 
(more blue centered) have lower log(G/S) than predicted. 
The strength of the correlation between residuals in 
log(G/S) and Ag-r is much lower than the other two pa¬ 
rameters (R = —0.15) with a larger probability of no cor¬ 
relation (P = 5x10“'’’). We test whether stellar mass is 
impacting this correlation between residuals in log(G/S) 
and Ag-r, by removing the correlation between residuals 
in log(G/S) and stellar mass. Running a Spearman Rank 
test on the correlation between the stellar mass corrected 
residuals in log(G/S) and Ag-r, we find that both the 
correlation strength becomes stronger (R = —0.25) and 
the probability of no correlation becomes several orders 
of magnitude smaller (P = 8x 10“^°). The correlation be- 
t ween An-r and stel lar mass has been shown previously 
in iStark et al.l (|2013fl and it is apparent that the stellar 
mass affects how this third parameter relates to residuals 
in log(G/S). The reliability of the measurement may be 
suspect for the smallest galaxies, since it requires mea¬ 
suring colors within i? 5 o,r- and between i? 5 o,r and Rr 5 ,r- 

The quantity that best meets our three criteria is b/a 
since it exhibits a strong correlation with the residuals 
in log(G/S) that is not affected by the correlation with 
stellar mass. The measurement of b/a is also the least 
likely to be affected by systematics from the photometry, 
including the convolution of galaxy SDSS and UKIDSS 
images to a common psf since it is a measure of the axial 
ratio in the outer disk. 

5.2. Physical Drivers of Residual Correlations 

For the purpose of defining PGF calibrations, we do 
not strictly need to understand the drivers of the correla¬ 
tions, but we explore them briefly here, deferring a more 
thorough discussion to future work. To aid our interpre¬ 
tation, we analyze log(G/S) residuals as a function of the 
best third parameter options and galaxy morphology in 
Figure [TOl Morphologies come from visual classification 
for the RESOLVE data set (Moffett et ah, submitted, 
Kannappan et al. in prep.). Since galaxy morphology 
is generally not available for large data sets and can be 
subjective and unreliable especially for small and edge- 
on galaxies, we have not used morphology for the general 
analysis of third parameters for improved gas mass esti- 
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mation. In Figure [TU] we plot galaxy morphology vs. 
Mr. 50 : and Ag_r and color code the residuals in 

log(G/S) such that green represents galaxies for which 
the PGF calibration underpredicts values of log(G/S), 
yellow represents galaxies for which the PGF calibration 
predicts values of log(G/S) similar to the measured val¬ 
ues, and red represents galaxies for which the PGF cali¬ 
bration overpredicts values of log(G/S). 

First we examine the relationship between morpho¬ 
logical type and /ir.BO in Figure [TUh . It appears that 
among late-type galaxies (T>0), the PGF calibration 
overpredicts values of log(G/S) for those galaxies with 
/ir .50 brighter than 22 and underpredicts log(G/S) for 
those galaxies with /ir, 5 o fainter than 22. For early- 
type galaxies, which mostly have fir ,50 brighter than 
22, we generally overpredict log(G/S). The transition 
around fir ,50 ~ 22 seems to indicate that even for a spe¬ 
cific type of galaxy, those galaxies that are lower sur¬ 
face brightness are on average more gas-rich than their 
higher surface brightness counterparts. Gas richness has 
long been associated wit h lo w surface brightness galaxies 
le.g. lde Blok et alllI996l and lO’Neil et al. 1200411 . and we 
show this result for a statistical population of galaxies 
using RESOLVE-A. 

Next we examine the relationship between morpholog¬ 
ical type and b/a in Figure [TUb . Within the population 
of late-type spirals (0 < T < 8), edge-on spirals tend 
to have underpredicted values of log(G/S). This phe¬ 
nomenon makes sense for larger edge-on spirals, which 
are observed to be redd er than their intrin sic colors due 
to dust extinction fe.g.. lGonrov et al1l20ICll show that for 
star forming galaxies over the narrow stellar mass range 
9.5 < log Mstar < 10, those with 6/a ~ 0.35 are red¬ 
der than those with b/a ^ 0.95). Thus the underpre¬ 
diction may result because these edge-on late type spi¬ 
rals are shifted redward of their true color measurement, 
where the typical value of log(G/S) is smaller^ The 
trend towards underpredicting edge-on galaxy values is 
less noticeable for edge-on irregular dwarf-type galaxies 
(T > 8) that are likely pure disk, low surface brightness 
galaxies with little int ernal extinction due to dust (e.g., 
iDalcanton et al.ll2004f l. For early-type galaxies, we see 
that in general values of log(G/S) are overpredicted with 
no obvious trend with b/a. 

We show the PGF calibrations using both model 
(green) and de-extincted (purple) u — J colors in Fig¬ 
ure EJi (we also show the PGF calibrations using g — r 
color in Figure ITTbl. The relationship between log(G/S) 
and color becomes steeper (less predictive) when using 
the de-extincted rather than m odel colors (FigurefTTh /cl. 
similar to the result shown in iJaskot et al.l ( 20I5ll . We 
do not find, however, that that the PGF correlation com¬ 
pletely disappears when using de-extincted colors, sug¬ 
gesting that while dust plays a role in the PGF calibra¬ 
tion, long-term star formation is also important (KI3). 
We suspect that star formation history is in fact related 
to dust extinction, which we plan to follow up in future 
work. In Figure ITTb we show the residuals in log(G/S) 
from the de-extincted PGF correlation as a function of 

® The reader may be concerned that a bias in stellar mass esti¬ 
mation as a function of extinction could lead to a spurious trend 
with axial ratio, but we have shown t hat our stellar masses are 
unbiased with respect to axial ratio in 


axial ratio, and we find that the strength and significance 
of the third parameter correlation are smaller but still 
highly significant. Additional effects, perhaps related to 
the morphology-surface brightness correlation, may also 
be important. 

Nevertheless, as the goal of this work is to provide em¬ 
pirical PGF calibrations for predicting gas masses, we 
prefer not to correct for extinction, since uncorrected 
colors provide the best (most linear and most predic¬ 
tive) calibrations. Moreover, uncorrected colors are not 
as dependent on modeling. 

Lastly we examine the relationship between morpho¬ 
logical type and stellar mass corrected Ag_r in Figure 
[IQI:. We correct Ag-r values for the corr e lation between 
Ag-r and stellar mass as in iStark et all (|20I3f) to show 
the effects more strongly. Late-type galaxies are gener¬ 
ally more red-centered with a slight trend towards bluer 
centers for later types (larger values of T), while early- 
type galaxies tend to be more blue centered as most o f 
these are blue-sequence E/SOs (|KannaDDan et all l2009tl 
because normal red-sequence E/SOs are excluded by our 
color cut. 

Within the late-type population, for a given morpho¬ 
logical type more blue centered galaxies have overpre¬ 
dicted values of log(G/S). This result is intriguing be¬ 
cause among late-type galaxies, blue centered galax¬ 
ies are associated with recent star forma tion events 
due to interactions (jKannapDan et aP 120041 1 . and have 
higher Mh.,/Mht ratios and lower Muj/Mstar ratios 
(jStark et al.ll20I3h . Including the molecular gas for such 
blue-centered late-type galaxies that are offset low in the 
color-limited PGF calibration may bring their total gas- 
to-stellar mass ratio (where total gas mass = I.4M^f/ 
+ M.H 2 ) in line with the general relationship between 
log(G/S) and color as shown in Figure 8 of KI3. 

The analysis of log(G/S) residuals as a function of 
galaxy morphology and the aforementioned three pho¬ 
tometric parameters reveals physical trends for all three 
that may have implications for galaxy evolution. The two 
parameters fir, 5 o and Ag_r, however, are both covariant 
with stellar mass, so their ability to reduce scatter is 
partially artificial. Using b/a allows us to minimize scat¬ 
ter in a physically meaningful way without introducing 
covariance into the PGF calibration. 

5.3. Modified Color-Limited PGF Calibration 

To use 6/a in PGF calibrations, we compute a linear 
combination of color and 6/a to be the new predictor. We 
call this combination of color and 6/a “modified color.” 
First, we fit a plane in color, 6/a, and log(G/S), min¬ 
imizing scatter in log(G/S) and only using those data 
points that fall within the blue color trim and red color 
cutoff for a particular color choice. The planar fit deter¬ 
mines coefficients such that modified color = rngX color 
-I- mix(6/a). 

We then use the relationship between modified color 
and log(G/S) to create tighter modified color-limited 
PGF calibrations. We use the same procedure as in (jH 
First we limit our data set to the 744 galaxies with reli¬ 
able HI data. Then we define a red modified color cutoff 
and a blue modified color trim just as in (21 Lastly we 
perform an ordinary least squares forward fit. Again we 
choose not to weight the fit by the errors on log(G/S) 
due to correlations between the uncertainties and both 
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Figure 10. Residuals in lQg(G/S) (dot col or) shown within plots of galaxy morphological type vs. (a) /i.r,50? (t>) and (c) stellar mass 
corrected Ag—r ffollowing IStark et al.ll2013l) . The color scale is shown at the top, where green marks galaxies for which the PGF calibration 
underpredicts values of log(G/S), yellow marks galaxies for which the PGF calibration values of log(G/S) similar to the measured values, 
and red marks galaxies for which the PGF calibration overpredicts values of log(G/S). Late-type galaxies (T > 0) with underpredicted 
values of log(G/S) tend to be lower surface brightness, more edge-on, and have redder centers (smaller values of stellar mass corrected 
Ag—r)’ These trends become somewhat less noticeable for dwarf types (T > 8), and in general for early-type galaxies values of log(G/S) 
tend to be overpredicted (with the caveat th at mo st red-sequence early-types are redder than our red color cut so not shown). We discuss 
the physical significance of these patterns in ESI 
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Figure 11. Test of the effect of dust on the relationships between log(G/S), color, and axial ratio, a) Log(G/S) vs. u — J color comparing 
the model (green) vs. de-extincted (purple) colors. Contours are shown encompassing 90% of all galaxies excluding confused galaxies and 
weak upper limits. Dashed lines show the color-limited linear fits to the data. The relationship steepens when using the de-extincted colors, 
b) Versus 6/a, the residuals in log(G/S) from the de-extincted u — J PGF calibration have a smaller strength of correlation and significance 
than the residuals computed from the model u — J PGF calibratio n, but the correlati on is still highly significant (t he red-dashed line shows 
the linear fit), c) Same as panel a except using g — r color (as in IJaskot et al.ll2Q15l') . Unlike Daskot et al.l l|2Q15ll we do not find that the 
relationship between log(G/S) and g — r disappears completely when using de-extincted colors, although we do see that the relationship 
steepens, meaning dust is a contributing factor to the relationship. 


galaxy color and log(G/S) itself. Table [3] gives the mod¬ 
ified color coefficients ttiq and mi, the slope of the linear 
fit, the offset in log(G/S), and the a of the modihed color 
PGF calibrations. The blue modified color trim and red 
modified color cutoff are also given as well as the number 
of galaxies used in each modified color fit. We note that 
the scatter for the PGF calibrations using modified col¬ 
ors based on u and g are typically reduced by ~0.03 dex 
as compared to the PGF calibrations using color only. 
Another advantage of using the modihed color calibra¬ 
tions is the generally increased baseline (or range) of the 
predictor value. 

We note that many other works have used planes 
dehned by color and other galaxy structural param- 


eters 

to define tighter PGF calibrations (IZhang et al. 

2009 

ilCatinella et al.ll2010HLi et al.ll2012HGatinella et al. 

2013 

), often choosing stellar mass surface density or sur- 


face brightness, which reduc e scatter partial l y due to co- 
variance with log(G/S). In iCatinella et al.l (j2013f l. the 
authors use the term “gas fraction planes” to refer to 
these PGF calibrations. We choose to use the terminol¬ 
ogy “modihed color” to refer specihcally to the predict¬ 
ing quantity, which is a linear combination of color and 
a structural parameter, in this work chosen to be b/a. 

6. PROBABILITY DENSITY FIELD PGF CALIBRATIONS 

Due to the limitations of the color-limited linear hts, 
we present in this section a new PGF calibration method 
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Table 3 

Modified Color-limited Photometric Gas Fraction Calibrations 


modified color 
(mag) 

slope 

(dex/mag) 

log(G/S) offset 
(dex) 

a 

(dex) 

blue trim 
(mag) 

red cutoff 
(mag) 

N galaxies 

1.822(m 

- r)”" -1- 0.609(b/a) 

-0.958 

3.038 

0.296 

2.0 

4.1 

580 

1.454(« 

-i)™- + 0.595ib/a) 

-0.989 

2.887 

0.291 

1.8 

3.8 

579 

1.140(u 

- J)"* -b 0.594(6/a) 

-0.981 

3.659 

0.304 

2.5 

4.5 

581 

1.075(u 

- bs')’" - 1 - 0.605(6/a) 

-0.952 

4.228 

0.302 

3.3 

5.0 

525 

3.563(5 

- r)"* + 0.534(t)/a) 

-1.002 

1.813 

0.281 

0.7 

2.6 

582 

2 . 444(5 

- i)"* + 0.550(6^) 

-0.984 

1.881 

0.287 

0.8 

2.6 

568 

1 . 592(5 

- J)™- -b 0.550(6/a) 

-0.980 

3.218 

0.313 

2.1 

3.8 

540 

1 . 435(5 

- K)^ + 0.618(6/a) 

-0.864 

3.678 

0.306 

2.8 

4.5 

488 


that improves over the color-limited fits in three ways. 
First, the new method provides a full probability dis¬ 
tribution in log(G/S) for each galaxy given its color or 
modified color rather than a single number for log(G/S). 
Second, we are able to model the galaxy distribution past 
the red color cutoff region where quenched galaxies co¬ 
exist with star-forming galaxies. Third, in the spirit of 
survival analysis, we develop a Monte Garlo method to 
reinsert the 22% of RESOLVE-A galaxies with unreli¬ 
able HI data in the calibration, thus determining the 
calibration for the entire volume-limited data set. We 
describe the full probability distribution method in HB.ll 
and we outline the Monte Carlo method to reinsert un¬ 
reliable galaxies in 116.21 Instructions on how to use the 
calibration are provided in 116.31 Lastly in 116.41 we exam¬ 
ine whether the model changes significantly if we add in 
the small subset of galaxies that have yir^tot fainter than 
— 17.33, but still have Mbary > 10®-^ Mq. For simplicity, 
throughout these sections we use the variable “me” to 
refer to either color alone or modified color that includes 
a h/a term. 


6.1. Calibration Description 

To start we use the same RESOLVE-A galaxies as in 
0 those with reliable detections or strong upper limits. 
Detections are considered reliable if they are not con¬ 
fused or are deconfused with systematic errors < 25% of 
the HI mass measurement and have S/N > 5. Limits 
are considered strong if they provide a gas mass mea¬ 
surement of Mgos < O.OSMstar- We show log(G/S) vs. 
modified {u — J)™ color for the RESOLVE-A data set 
limited to reliable gas data in Figure [12^. This defini¬ 
tion excludes 211 galaxies (of 955 total) in RESOLVE-A, 
however we describe how a technique based on survival 
analysis principles to include these galaxies in the PGF 
calibration in > 16.21 

For galaxies with upper limits, we cannot know the 
value of the HI mass, only that it is less than a certain 
value. We choose to set all galaxies, both detections 
and limits, with log(G/S) < —1.3 to log(G/S) = —1.3, 
which is 5% of the stellar mass in linear units. This 
choice is fine for measuring a galaxy’s baryonic mass (as 
in the companion paper Eckert et al. in prep.), because a 
galaxy’s baryonic mass will not be significantly affected 
if it is at most 1.05Mstaf In Figure [T2h. we demonstrate 
the data replacement with red arrows. 

Next we determine the density field of galaxies in 
log(G/S) vs. me by binning in both dimensions and mea¬ 
suring the number of galaxies in each 2D bin or cell of 
log(G/S) and me. The density field is shown in Figure 


[Hd. We use bin sizes of 0.2 dex in log(G/S) and 0.2 mag 
in (u — J)"*. When determining the bin size for each me, 
we start with bin sizes of 0.2 mag and require a minimum 
of 10 bins in color to ensure adequate sampling of the cal¬ 
ibration. If there are less than 10 bins (as is the case for 
colors with small ranges such as g — r) we switch to bin 
sizes of 0.1 mag. 

To model the PGF probability density field we assume 
two distinct populations: detections and upper limits. 
The detections lie on a line with Gaussian scatter that 
widens towards redder colors. The limits appear only at 
red me and are confined to a value of log(G/S) = —1.3 
although they may in fact have much lower values. The 
model consists of nine parameters given by Aq to Ag. 

To model the population of detections we first model 
the peak value of the Gaussian at each me, designated 
as po, as a log-normal function of me as given in equa¬ 
tion [TJ Aq provides the normalization, Ai is the location 
parameter in me, and A 2 is the shape parameter. 


Po 


(mcA^V^) 2A2 ^ 


( 1 ) 


The number density of detections {D) is then described 
by a Gaussian with peak value po that changes with me, 
with mean log(G/S) value {A^mc -b A 4 ) that decreases 
linearly for larger values of me, and with standard devia¬ 
tion A^mc that widens towards larger values of me. The 
model for the detections is given in equation [2j 


D = Po X exp( 


-{log{G/S) - (Apme + A^))'^ 
2(A5mc)2 


( 2 ) 


At red me, there are galaxies with values of 
log(G/S) < —1.3 that belong to the detection rather 
than limit population. Because we have chosen to set 
all galaxies with values of log(G/S) < —1.3 to be equal 
to —1.3, we must divide up equation [5] to account for the 
tail of the detection population with log(G/S) < —1.3. 
For values of log(G/S) > —1.3, the model is as given by 
equation[2J but for values of log(G/S) = —1.3, we model 
the detection population by integrating the tail of the 
Gaussian at a given me over all log(G/S) < —1.3. Thus 
the detection population for log(G/S) < —1.3 is given in 
equation [3l 



er/( 


[-1.3- (Asme + Aj)] 

\/2(A5mc)2 

( 3 ) 
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Figure 12. Log(G/S) vs. modified {u — J)™' color for RESO LVE- A. The modified color is the linear combination of {u — J)"* and b/a 
that produces the best plane fit with log(G/S) as described in 115.31 a) Scatter plot showing reliable HI detections with small black dots, 
and strong upper limits with downward black arrows. Galaxies with log(G/S) values < —1.3 (marked by the black line), whether they are 
detections or limits, are replaced with —1.3 (as indicated by the red arrows), b) Contour plot showing the PGF density field of log(G/S) 
vs. modified {u — j)”* color, where each cell’s value is the number of galaxies in that cell. The contou rs ar e spaced finely over 100 intervals 
between 0 and 34 (the maximum number of galaxies found in a cell). We fit the model described in llfi.ll to this PGF density field. 


The upper limit population is evaluated only at 
log(G/S) = —1.3 and is modeled as a Gaussian function 
of me with a peak value Ag, color shift Ay, and standard 
deviation Ag as given in equation |4l The Gaussian shape 
captures the fact that blue galaxies rarely have low gas 
fractions and on the other side that there is a natural 
decline in galaxies towards redder colors. 

L = AeX exp( (4) 

The full-probability PGF model is the combination 
of equations El El and m and we fit this model to 
the PGF density field using the MPFIT2DFUN pack¬ 
age (jMarkwardtl EPOOTFI . whi ch performs a Levenberg- 
Marquardt least squares fit (iMord I1978D . We weight 
each cell value in the density field, which is the num¬ 
ber of galaxies in the cell N, by 1/N. From the fit we 
obtain the parameters Aq through Ag that best describe 
the PGF calibration. We recognize that the probabil¬ 
ity density field model is more complex than a simple 
linear fit. However, the advantages of being able to fit 
both the detection and limit populations and of being 
able to reinsert the data from galaxies without adequate 
HI data (see 96.21) make it a more powerful technique for 
estimating gas masses (see ©• 

6.2. Reinserting the Missing Galaxies 

To measure probability density field PGF calibrations 
with the complete RESOLVE-A data set, we implement 
a custom version of survival analysis following its ba¬ 
sic algorithm: create a model of the missing/censored 
data points based closely on the available data and sam¬ 
ple it with Monte Garlo methods. The model for each 
censored/missing data point is determined from the un¬ 
censored data in an initial fitting round as in traditional 

® http://purl.com/net/mpfit 


survival analysis, with the restriction that limits must be 
drawn from a model distribution truncated at the limit 
values. We then iterate this procedure a second time 
using the updated model resulting from the first round. 
The steps are as follows. 

1) After performing the initial PGF model fit from 96.11 
using only the reliable HI data in RESOLVE-A, we ob¬ 
tain a PGF model with initial best fit parameters A^^q 
through Ai^g. Using these initial best fit parameters, 
we obtain a probability distribution in log(G/S) for each 
galaxy with unreliable HI data given its me. For galaxies 
with weak upper limits, we restrict the probability distri¬ 
butions in log(G/S) to only those values below the upper 
limit, then we renormalize the probability distribution. 
For low S/N detections and confused galaxies we do not 
place any restrictions on the probability distribution in 
log(G/S). 

2) We then randomly assign a value of log(G/S) for 
each of these 211 galaxies from each respective proba¬ 
bility distribution using the inverse transform sampling 
method. 

3) Next, we create a new PGF density field that in¬ 
cludes both the 744 data points with reliable HI and 
the 211 data points with randomly assigned values of 
log(G/S) from their respective probability distributions. 
We fit the PGF calibration model as described in gnu 
and save the best fit parameters A„^o through A„_g. We 
also compute residuals for each best fit model as well as 
the reduced goodness of fit, assuming the variance on 
the data value for each cell is cr^ = N (the number of 
galaxies in the cell). 

4) We perform steps 2 & 3 100 times and save the 
best fit parameters and reduced values for each round. 
We then calculate the median values of the parameters 
based on those 100 rounds Amed,o through Amed,8, which 
describe the new model for the PGF density field. We 
save the median reduced as well. 
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5) We iterate once more through steps 1-4 using the 
model parameters through Amed,8 to produce the 

initial probability distributions in log(G/S) for the 211 
galaxies with unreliable HI data. This iteration allows 
us to determine log(G/S) probability distributions based 
on the PGF density field calibrated with all RESOLVE-A 
galaxies, not just those with reliable HI data. After this 
second iteration, we take the values of Amed.o through 
Amed.s to be the final best fit parameters Af o through 
A/8 for the full-probability PGF model. We have tested 
whether more iterations are necessary, however after the 
second iteration the values of Amed.o through A^ed.s 
change very little (< 3% of the final value). The final 
values do change significantly from the initial fit (which 
includes only the reliable HI galaxies), with final values 
up to ^20% different from the initial values. 

Figure [13] shows a contour plot of the best fit model 
and residuals using modified {u — J)'” color, where the 
residuals have been normalized by the number of galax¬ 
ies in each cell. The 211 galaxies with unreliable HI data 
have been assigned values of log(G/S) pulled at random 
from their respective distributions in log(G/S), as in step 
2, using the final best fit parameters. The residuals are 
most significant for the highest log(G/S) cells at red mod¬ 
ified color where there are few galaxies (often just one per 
cell). 

6.3. How to Use the PGF Calibration 

The parameters in Tables ID and [5] provide the best 
fit models to the PGF density field probability distribu¬ 
tion, which is not normalized at a given value of me. To 
create conditional probability distributions in log(G/S) 
at each me, the PGF distribution for a given me from 
equations laisi and m should be divided by the integral 
of the PGF at that mcQ Examples of such distributions 
are shown in Figure [T4l for modified {u — J)™ color me 
= 3.1 (blue), 4.4 (green), and 5.3 (red). We note that 
the sharp spike in the log(G/S) probability distribution 
is entirely due to the way in which upper limits and low 
gas-fraction galaxies are treated in the model and should 
not be interpreted literally as probability distributions 
for this population. 

6.4. Adding in High M^ary/L Galaxies 

As a check, we have also performed the full-probability 
PGF calibration using a slightly altered RESOLVE-A 
data set. This data set aligns with the spectroscopic 
observing data set for the RESOLVE survey. We con¬ 
sider all galaxies with Mr^tot brighter than —17.33 as well 
as any galaxies with M^ary > 10®'^ Mq no matter their 
thus including the highest baryonic mass-to-light 
ratios that fall out of the absolute magnitude limited 
sample. 

For this test we use the PGF calibrations from TablesjH 
and [5] to estimate gas masses for the galaxies with Mr^tot 
fainter than —17.33 that do not have adequate gas data. 
Using the median value of log(G/S) from each proba¬ 
bility distribution we estimate the gas mass for these 
galaxies. We then compute M^ary for each galaxy in the 

^ We provide an IDL code at 

https://github.com/keckert7/codes/pred_loggs_dist.pro that 

looks up the best fit parameters for a given me and provide the 
probability distributions for any set of galaxies with provided me. 


RESOLVE-A region, adding the stellar mass to either the 
measured M^as for those galaxies with reliable detections 
and strong upper limits, or the predicted for those 
galaxies with low S/N detections, confused profiles, or 
weak upper limits. Adding galaxies with Mr^tot fainter 
than -17.33 and Mbary > 10® ® Mq yields a RESOLVE-A 
data set with 39-40 additional galaxies depending on the 
choice of PGF calibration. Thirty-nine of those galaxies 
have reliable HI data. 

With these galaxies added, we c an no w ru n th rough 
the same procedure as detailed in 96.11 and 96.21 The 
parameters for the best fit models are given in Tables 
in and |7| The PGF calibrations yield similar log(G/S) 
distributions with or without the high baryonic mass-to- 
light ratio galaxies. In Figure |T4| grey solid lines show 
the log(G/S) probability distributions for me = 3.1, 4.4, 
and 5.2 when including the high baryonic mass-to-light 
ratio galaxies. The largest difference is seen for the bluest 
galaxies, where the high baryonic-mass-to-light galaxies 
very slightly shift the probability distribution towards 
higher values of log(G/S). 

7. DISCUSSION 

In this section we perform two comparisons of PGF 
calibrations from this work and other work using the 
RESOLVE-B data set. First, we compare predicted gas 
masses from PGF calibrations with measured gas masses 
from the RESOLVE-B data set as a function of stellar 
mass. To test the different PGF calibrations, we select 
only galaxies in RESOLVE-B that have reliable HI de¬ 
tections (294 out of 487), as upper limits are not ideal 
candidates for this test. Second, we compare the actual 
RESOLVE-B distribution of log(G/S) in bins of {u - J)"* 
color with the predicted distributions of log(G/S) based 
on these different calibrations. 

To test the probability density field model we use the 
modified (u — J)™ color calibration (top row of Figure 
HSj). Since the calibration provides a probability distri¬ 
bution of log(G/S) for each galaxy, we use the median 
value from the probability distribution in log(G/S) for 
the purposes of the first test. In the top panel of Figure 
US! we find that using the probability density field model 
PGF calibration yields a negligible offset and small scat¬ 
ter = 0.343 dex. For the lowest stellar masses, we see 
asymmetric scatter towards larger measured than pre¬ 
dicted gas masses. As indicated by the arrow in the 
bottom panel of Figure Hsj this asymmetric scatter is 
consistent with covariance between stellar mass, plotted 
on the x-axis, and the predicted gas mass in the denom¬ 
inator of the y-axis (note that estimated G/S is multi¬ 
plied by stellar mass to obtain gas masses). We also show 
the distribution of log(G/S) predicted for RESOLVE-B 
galaxies in four bins of (u — J)™ color, drawing a random 
value from the log(G/S) distribution of each RESOLVE- 
B galaxy. The probability density field model works well 
across all colors. 

To test our color-limited linear fits, we use the cali¬ 
bration based on {u — J)'" color. While the {u — J)"* 
calibration has higher scatter than calibrations based on 
optical colors, the {u — J)'" calibration has a larger pre¬ 
dictive range as compared to typical errors on color (con¬ 
sidering the slopes that would multiply the error; see Ta¬ 
bles |2| and |3| In the second row of Figure [151 we show 
the log of the measured gas mass divided by the pre- 
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Table 4 

Full-Probability Photometric Gas Fractions Calibrations for Color Only 


color 

(mag) 

^f,o 

^/,i 

7l/,2 

^/,3 

^/,4 

^/,5 



00 

reduced 

{u - r)”^ 

17.22 

0.33 

0.22 

-1.67 

2.56 

0.21 

47.65 

2.21 

0.14 

1.353 

(u - i)"* 

35.05 

0.46 

0.24 

-1.37 

2.40 

0.19 

68.44 

2.55 

0.20 

1.518 

(u - J)^ 

32.55 

1.00 

0.17 

-1.04 

3.07 

0.11 

53.39 

3.88 

0.26 

1.257 

(u — K)'^ 

30.18 

1.26 

0.13 

-0.99 

3.75 

0.09 

42.35 

4.79 

0.33 

1.422 

[g - rT 

19.70 

-1.15 

0.48 

-3.65 

1.51 

0.87 

115.1 

0.72 

0.05 

1.831 

(g - i)- 

19.99 

-0.70 

0.45 

-2.47 

1.58 

0.55 

72.29 

1.05 

0.09 

1.344 

(9 - jr 

16.99 

0.51 

0.19 

-1.52 

2.80 

0.18 

42.87 

2.37 

0.16 

1.036 

{g - 

31.08 

0.91 

0.15 

-1.22 

3.27 

0.14 

61.95 

3.29 

0.22 

1.246 


Table 5 

Full-Probability Photometric Gas Fractions Calibrations for Modified Color 

modified color 
(mag) 

^/,o 

-4/,i 

^/,2 

^f,3 

7l/,4 

^/,5 

^/,6 

^/,7 

^f,s 

reduced 

1.822(u - r)”* + 0.609(b/a) 

38.78 

1.05 

0.20 

-0.89 

2.80 

0.09 

47.95 

4.44 

0.29 

1.125 

1.454(m - i)"* -1- 0.595(b^) 

37.78 

0.97 

0.21 

-0.93 

2.70 

0.10 

51.88 

4.10 

0.27 

1.104 

1.140(m - J)"* -1- 0.594(b/a) 

38.13 

1.25 

0.16 

-0.93 

3.49 

0.08 

44.57 

4.82 

0.31 

1.066 

1.075(m - K)^ + 0.605(b/a) 

34.54 

1.43 

0.14 

-0.88 

3.92 

0.07 

40.53 

5.53 

0.34 

1.157 

3.563(g - r)"* -I- 0.534(b/a) 

42.65 

0.37 

0.37 

-0.97 

1.74 

0.18 

62.12 

2.92 

0.21 

1.166 

2.444(g - i)™- -|- 0.550(b^) 

40.95 

0.45 

0.36 

-0.99 

1.87 

0.17 

60.79 

2.95 

0.21 

1.228 

1.592(g - J)"* -1- 0.550(b/a) 

36.40 

1.10 

0.18 

-0.94 

3.08 

0.10 

51.50 

4.16 

0.26 

1.028 

1.435(g - it')’" -1- 0.618(b/a) 

33.79 

1.37 

0.14 

-0.91 

3.83 

0.08 

45.48 

5.10 

0.31 

1.182 


Table 6 

Full-Probability Photometric Gas Fractions Calibrations for Color Only Including High 

Mbary/L Galaxies 

color 

(mag) 

^f,0 


^/,2 

^f,3 

7l/,4 

^/,5 

^/,6 

71/.7 

^/,8 

reduced x^ 

(« - r-)’" 

34.56 

0.31 

0.23 

-1.66 

2.57 

0.24 

90.35 

2.21 

0.15 

1.514 

(u - i)^ 

35.01 

0.44 

0.24 

-1.44 

2.54 

0.20 

68.48 

2.55 

0.20 

1.474 

(u - J)^ 

32.37 

1.00 

0.17 

-1.10 

3.28 

0.12 

53.83 

3.88 

0.26 

1.308 

(u - K)^ 

29.47 

1.26 

0.13 

-1.03 

3.93 

0.10 

42.53 

4.78 

0.33 

1.446 

{g - rT 

19.24 

-1.17 

0.49 

-3.74 

1.56 

0.93 

115.8 

0.72 

0.05 

2.033 

{g - ir 

19.75 

-0.71 

0.45 

-2.54 

1.63 

0.58 

72.31 

1.05 

0.09 

1.514 

{g - JT 

32.44 

0.50 

0.20 

-1.54 

2.87 

0.22 

82.06 

2.39 

0.17 

1.098 

(g - 

29.46 

0.91 

0.15 

-1.29 

3.46 

0.15 

62.12 

3.29 

0.23 

1.323 


Table 7 

Full-Probability Photometric Gas Fractions Calibrations for Modified Color Including High Galaxies 


modified color 
(mag) 

^f,0 

7^/,l 

to 

^f,3 

^fA 

^f,5 

^f,6 

^f,7 


reduced x^ 

1.822(u 

- r)’" + 0.609(b/a) 

38.82 

1.05 

0.21 

-0.90 

2.84 

0.10 

48.06 

4.44 

0.29 

1.194 

1.454(m 

- O’" -1- 0.595(b^) 

37.44 

0.96 

0.21 

-0.96 

2.79 

0.11 

51.32 

4.10 

0.27 

1.188 

1.140(m 

- J)"* -1- 0.594(b/a) 

36.98 

1.24 

0.16 

-0.98 

3.66 

0.08 

44.90 

4.82 

0.31 

1.154 

1.075(m - 

- is:)’" -1- 0.605(b/a) 

33.11 

1.43 

0.14 

-0.91 

4.05 

0.07 

40.99 

5.53 

0.34 

1.219 

3.563(g 

- r)"* -1- 0.534(b/a) 

41.40 

0.36 

0.37 

-1.01 

1.82 

0.19 

61.93 

2.91 

0.21 

1.299 

2.444(g 

- i)™- + 0.550(b^) 

39.53 

0.44 

0.35 

-1.03 

1.97 

0.18 

60.35 

2.95 

0.21 

1.309 

1.592(g 

- J)’" -I- 0.550(b/a) 

35.19 

1.10 

0.18 

-0.97 

3.20 

0.10 

51.74 

4.16 

0.26 

1.171 

1.435(3 - 

- is')’" - 1 - 0.618(b/a) 

32.02 

1.37 

0.14 

-0.94 

3.96 

0.08 

45.97 

5.10 

0.31 

1.216 
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234523456 
1.140(u-J)'^ + 0.594(b/a) 1.140(u-J)^ + 0.594(b/a) 


Figure 13^ Best fit model to the PGF density field using modified (u — J)^ color, a) The best fit model is shown as a contour plot (as 
for Figure fT^ divided into 100 levels between 0 and 44, the maximum number of galaxies found in a cell after reinsertion of galaxies with 
inadequate HI data), b) Residuals for each cell normalized by the number of galaxies in each cell so the range is from —1 to 1. Overall the 
model performs well, although there are large residuals especially around the edges of the calibration where the number density of galaxies 
is low. 
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log(G/S) 

Figure 14. Probability distributions for three galaxies with mod¬ 
ified u — = 3.1, 4.4, and 5.2, which are labeled on the 

plot. For blue modified colors the distribution in log(G/S) is Gaus¬ 
sian (me = 3.1). The a value of the Gaussian widens for redder 
modified colors (me = 4.4), where the population of galaxies with 
log(G/S) < —1.3 also increases. The sharp spike at log(G/S) = 
— 1.3 is due to the way in which galaxies with upper limits or low 
gas-fractions are treated in the model. For the reddest colors (me = 
5.2)the contribution from the upper limit population with log(G/S) 
< —1.3 becomes dominant. In solid grey we show the probability 
distributions for the same modified colors when using the fits that 
include high baryonic mass-to-light galaxies as explained in 16.41 
The largest difference occurs for the bluest galaxy, although the 
difference is quite small. 


dieted gas mass as a function of stellar mass. We find 
a negligible offset between the measured and predicted 
gas mass measurements and a scatter of cr = 0.361 dex, 
slightly higher than found in the color-limited PGF cal¬ 
ibration itself and slightly higher than the probability 
density field model. We next show the log(G/S) distri¬ 
butions predicted for RESOLVE-B with the color-limited 
linear ht, where we add random scatter with cr = 0.31 
dex to log(G/S) estimates. The second test reveals the 
weakness of the color-limited approach as it cannot re¬ 
produce the distribution of log(G/S) for red galaxies in 
RESOLVE-B. 

In the third row, we test the calibration based on the 
NFGS given in K13, log(G/S) = -0.98 (m- J)”" -f 2.70. 
The NFGS has been run through the same photometric 
pipeline and stellar population modeling code as for this 
work, so we are able to use the same SED modeled colors 
and stellar mass estimates to predict gas masses. There 
is a signiheant offset (0.172 dex) towards underpredict¬ 
ing gas mass measurements, with that offset increasing 
for low stellar mass galaxies. The same issue can be seen 
in the second test, where the K13 calibration systemati¬ 
cally underestimates the distribution of log(G/S) for blue 
galaxies and overestimates the distribution of log(G/S) 
for red galaxies in RESOLVE-B. 

A similar problem is seen for the color- 
limited (NUV—r < 4. 5) PGF c ali bratio n 
based on GASS from iGatinella et al.l (j2013l l. 
log(G/S) =-0.2341og(/i,)-0.342(NUV-r) -k 2.329. 

This color-limited calibration attempts to reduce resid¬ 
uals for gas-rich galaxies. We note that the GASS 
survey calibration does not take into account a helium 


correction factor so we multiply the gas mass value 
obtained from this calibration by 1.4 to compare with 
the gas mass measurements in RESOLVE-B. To match 
their parameters we use NUV—r total colors measured 
for the RESOLVE survey, which are similar to the Kron 
NUV—r colors measured for the GASS survey, although 
subject to differences in the optical sky background 
subtraction and magnitude extrapolation. GASS defines 
= M»/(27ri?|Q , 2 ), where z-band Petrosian 

half-light radius in kpc. We use the stellar mass mea¬ 
surements from this work to nonetheless compute Mstar 
and convert predicted G/S to predicted gas mass0 Due 
to the constraint of requiring NUV total magnitudes 
we can use only 281 galaxies for testing the GASS 
calibration. The bottom row of Figure ITSl shows that as 
for the NFGS, there is a large offset (0.235 dex) between 
the predicted and measured gas masses that increases 
for lower stellar mass galaxies, in the sense that the 
calibration underpredicts gas masses. Also apparent is 
the underestimate of the distribution of log(G/S) for 
blue galaxies in RESOLVE-B. 

Both the NFGS and GASS are representative galaxy 
samples. The NFGS is a representative sample of galax¬ 
ies designed to cover all morphological types over a range 
of galaxy magnitudes in rough proportion to the R-band 
luminosity function (I Jansen et al.l 1200011 . The GASS 
sample is designed to examine the full range of massive 
gala xies with a flat distri bution in stellar mass > 10^° 
Mq (|Gatinella et al.l[201fl[ l. Both samples have complete 
HI data to fractional mass limits < 1-5% of the stel¬ 
lar mass, making them seem like ideal calibration sam¬ 
ples. They do not, however, statistically represent the 
galaxy population of the local universe, in particular the 
d ominance of low -mass gas-rich dwarfs. We note that 
in iLi et al.l (j2012tl , the authors create a calibration using 
GASS that includes NUV—r color, stellar mass, stellar 
mass surface density, and g — i color gradient and more 
successfully r eplicates gas m a sses f or gas-rich galaxies 
than does th e IGatinella et aP (|2012ll calibration. We do 
not show the ILi et al.l " 201211 calibration due to the large 
quantity of covariant variables required, which are also 
difficult to calculate from the data that we have available. 

We hav e also tested (but do not show) the PGF calibra¬ 
tion from iHuang et al.l (|2012ll which uses the ALFALFA 
survey to define a PGF calibration. This calibration re¬ 
lies on SDSS Petrosian g — r color and stellar mass sur¬ 
face density. The latter is defined in the same way as 
for GASS with z-band Petrosian radii. We hnd that the 
ALFALFA calibration overpredicts gas masses with an 
offset of —0.19 dex. This over prediction of gas masses 
is not surprising since the ALFALFA survey is HI flux 
limited and weighted towards gas-rich objects. 

Lastly, we have tested (but do no t sho w) the PGF 
calibr ations from iKannapoanl (|2004r i and I Zhang et al.l 
(j2009l l. both based on heterogeneous samples of HI detec- 
tions from the HyperLeda database (|Paturel et al.|[200^ 


® Th e GASS stell a r mas ses are computed using the SED fitting 
code of ISalim et al.l 1120071 ) that populates a g rid of stellar popu¬ 
lation models with a Chabrier IMF using the IBruzual Ik Chariot! 
{2003) stellar popul ation synthesis co d e. Comparison betwee n 
the stellar masses of jSalim et al.l II2005I1 , IKauffmann et ahl II2003I 1 , 
IKannappan et a P Il200 9ll. and K13 indicates that the stellar masses 
computed as in l^ahni^t al.l II2007I 1 should be extremely similar to 
the K13 stellar masses used in this work (within 0.1 dex). 
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crossmatched to c atalogs from SDSS and 2MASS. The 
iKannappanl (I2004D ca librations rely on c atalog u — r or 
u — K color, while the lZhang et all l)2009f) calibration re¬ 
lies on a combination of catalog g — r color and i-band 
surface brightness. Both use s tellar mass e s com puted 
following the prescription from iBell et (1200,1( 1 . We 

mimic these quantities, removing galaxies without qual¬ 
ity 2MASS AT-band data (limiting t he analysis to 182 
gala xies). ThePGF cal ibrations from iKannappanI (|2004D 
and iZhang et al.l (|2009H show very large scatter of ~0.7 
and ~0.63 respectively, but produce surprisingly small 
offsets (—0.06 and —0.02 dex respectively), likely due to 
fortuitously cancelling systematics (underrepresentation 
of the dwarf galaxy population but overrepresentation of 
galaxies with high 21cm flux). 

Based on these comparisons, we stress the importance 
of matching the selection criteria and completeness of the 
calibration data set to the properties of the data set for 
which predictions are desired. Using PGF calibrations 
for a small representative galaxy sample cannot repro¬ 
duce the gas masses for a volume-limited survey, which 
is dominated by low-mass, gas-rich galaxies. 

The purpose of these z=0 PGF calibrations is to enable 
measurement of galaxy cold gas masses for large surveys 
of the nearby universe. Previous uses of PGF calibrations 
include identifying gala xy mass scales con nected to tran¬ 
sitions in gas-content (IKannap panl 1200 4 K13), estima¬ 
tion of the HI mass function ((Zhang et al.ll2009H . explo¬ 
ration of the dependence of t he stellar mass-meta llicity 
relationship on gas content (jZhang et al.l I2009H . and 
galaxy clustering sta tistics within galaxy populations of 
varying gas content (|Li et al.ll2012h . 

These new calibrations are suited for similar studies 
in volume-limited nearby galaxy surveys, for which ob¬ 
taining complete HI data is not feasible. In a companion 
paper (Eckert et al. in prep.), we use the probability den¬ 
sity field model to create distributions of log(G/S) for 
RESOLVE and ECO galaxies, which we then combine 
with the stellar mass likelihood distributions produced 
by the SED fitting described in il1.2l to produce baryonic 
mass likelihood distributions for each galaxy. We use 
Monte Carlo sampling from these baryonic mass likeli¬ 
hood distributions to construct baryonic mass functions 
for the RESOLVE and ECO data sets. In Florez et al. 
(in prep.) the probability density field model is used to 
examine the gas content of galaxies in the lowest-density 
environments in ECO. Finally, we caution against the use 
of these calibrations at much higher redshift as galaxy 
properties (including gas-content and color) may exhibit 
different relationships at higher redshifts. 

8. CONCLUSIONS 

We have presented the newly reprocessed 
NUVugrizVJiJA" photometry for the RESOLVE 
survey, all provided in a machine readable table (Table 
[1]). We have also provided new z=0 volume-limited 
calibrations of the PGF technique using both linear 
fits (limited to exclude the reddest colors) and, more 
powerfully, a 2D model of the full probability density 
field distribution of log(G/S) vs. color (or “modified 
color,” including axial ratio). We highlight the main 
conclusions of this work below: 

• The new photometry uses improved background sub¬ 
traction and multiple flux extrapolation methods to yield 


brighter magnitudes, bluer colors, and larger sizes than 
catalog photometry (see ^3.11 and Figures |3] and 0]). 
Multiple flux extrapolation techniques allow us to mea¬ 
sure systematic errors, key for estimating reliable stellar 
masses (see T3.2I1 and computing SED-modeled colors. 

• This new photometry reveals a real increase in scatter 
around the red sequence, which we attribute to the fact 
that our flux extrapolation routines do not suppress color 
gradients, unlike the SDSS model magnitude algorithm 
(see Figure |4|). 

• We provide linear PGF calibrations between 
log(G/S) and color, however due to the breakdown in the 
relationship between log(G/S) and color for red galaxies, 
these linear fits must be color-limited and do not per¬ 
mit estimation of gas masses for galaxies beyond the red 
color cutoff (see 0and Table [2|). The color-limited PGF 
calibrations also do not account for the 24% of galaxies 
that are confused, have weak upper limits, or lack reliable 
detections, biasing the fits towards higher gas-to-stellar 
mass ratios at the red end. 

• We find that axial ratio correlates significantly with 
residuals in log(G/S), independent of stellar mass, which 
has a covariant correlation that is also artificially en¬ 
hanced for an approximately baryonic mass limited sam¬ 
ple such as ours (see >15.11 and Figures 151 and ITU)) . The 
correlation between log(G/S) residuals and axial ratio 
may be at least partially related to dust and internal 
extinction, which also may partially drive the PGF cal¬ 
ibration along with long term star formation rate. Pla¬ 
nar fits between color, axial ratio, and log(G/S) provide 
linear combinations of color and axial ratio (“modified 
color”) that yield tighter PGF calibrations than using 
color alone (see > 15.31 and Table [3|). 

• We provide a new type of PGF calibration using a 
2D model to fit the entire probability density field of 
log(G/S) vs. modified color (see >J6] and Figures [12] and 
ini). Within this calibration we are able to statistically 
model galaxies missing from the color-limited linear fits. 

• The probability density field model PGF calibra¬ 
tion yields log(G/S) probability distributions for indi¬ 
vidual galaxies (see Figure m. For red galaxies, for 
which linear fits cannot accurately predict log(G/S), the 
full-probability PGF yields two-component distributions, 
representing both quenched and star forming galaxies 
(see Figure ITSl). 

• We test our color-limited and full-probability PGF 
calibrations as well as literature calibrations using the 
RESOLVE-B data set and find that our calibrations are 
well suited to predicting gas masses for a volume-limited 
sample (see Figure 1151) . Previously published calibra¬ 
tions generally systematically over- or underpredict gas 
masses, partly due to the different selection criteria of 
the calibration samples. 

We emphasize the importance of defining a calibration 
sample with the same selection criteria and completeness 
as the sample for which predictions are desired. Here 
we have presented new PGF calibrations using the RE¬ 
SOLVE survey that are ideal for application to volume- 
limited data sets complete to a baryonic mass (or optical 
luminosity) limit. These PGF calibrations will be used to 
determine baryonic masses in a companion paper (Eckert 
et al. in prep.) that examines the baryonic mass function 
as a function of environment in RESOLVE-B and ECO . 
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Figure 15. Comparison of PGF calibrations from this work and the literature using 21cm data for RESOLVE-B. Each row shows a 
different calibration. First Column: Log of measured over predicted gas mass for RESOLVE-B galaxies with HI detections as a function 
of stellar mass. For the probability density field model (top row), we use the median value of each galaxy’s log(G/S) distribution as the 
single value gas mass prediction, and we find a negligible offset and and small scatter in the residuals. The color-limited linear fit from 
this work (2nd row) also p roduces a negligible offset, but slightly larger scatter in the residuals. The color-limited linear fits from K13 
and ICatinella et al.l II2013I) (3rd and 4th rows) both yield large offsets towards underpredicting gas masses. The asymmetric scatter for 
low stellar mass galaxies is caused by the covariance betwee n stellar mass and the estimated ga s ma ss in the denomina tor of the y-axis 
(arrow in bottom panel). Comparisons to other calibrations IIKannappanll2Q0^. IZhang et al.ll2009l. and IHuang et ani2012l) are discussed in 
the text (0. 2nd-5th Columns: Actual (green) and estimated (purple) distributions of log(G/S) for RESOLVE-B, where all values of 
log(G/S) < —1.3 are set equal to —1.3. For the probability density field (top row), we draw random values from each galaxy’s log(G/S) 
distribution. The estimated log(G/S) distributions are consistent with the data for all colors. For the color-limited linear fits, we add 
random scatter with cr ~ 0.3 to the log(G/S) estimates for each galaxy. The color-limit ed linear fit from this work (2nd row) performs well 
for blue galaxies, while the color-limited linear fits from K13 and fGatinella et al.l II2013I1 (3rd and 4th rows) underestimate the distributions 
of log(G/S) for blue galaxies. None of the color-limited linear fits performs well for red galaxies, necess itating the probability density field 
model. We attribute the underprediction of gas masses for low mass and blue galaxies by the K13 and IGatinella et al.l 1120131 ! calibrations 
to their samples (NFGS and GASS respectively) not being volume-limited samples of the galaxy population in the nearby universe. 
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